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SECTION  I 
INTRODUCTION 

The  response  of  aircraft  structures  to  the  blast  overpressure 
generated  by  a  nuclear  explosion  is  an  important  consideration  of  the 
nuclear  survivability  for  aircraft  during  several  phases  of  its  mission. 
The  computer  code  NOVA  (Nuclear  Overpressure  Vulnerability  Analysis) 
given  in  References  1  and  2  was  developed  by  Kaman  AviDyne  for  the  Air 
Force  Weapons  Laboratory  (AFWL)  originally  to  predict  the  dynamic 
response  of  individual  aircraft  structural  elements,  such  as  stringers, 
frames  and  unstiffened  panels,  to  the  transient  overpressure  loads 
associated  with  a  nuclear  burst.  The  N0VA-2S  version  (Reference  3) 
extended  the  panel  response  subroutine  DEPROP  (Dynamic  Elastic  Plastic 
Response  of  Panels)  to  handle  flat  or  circular  curved  stiffened  panels 
in  both  the  elastic  and  inelastic  response  regimes.  Thus,  N0VA-2S  could 
analyze  stringers  and  frames  in  combination  with  the  panel  skin  and, 
thereby,  represent  major  components  of  an  aircraft,  such  as  fuselage 
sections  and  stiffened  panels  between  spars  and  ribs  of  the 
wing,  vertical  tall  and  horizontal  stabilizer.  While  the  N0VA-2S 
computer  code  is  used  for  nuclear  overpressure  loading,  the  special 
version  N0VA-2LTS  Is  available  for  general  pressure  loadings  and,  in 
particult  ,  tor  use  with  measured  pressure  test  data  stored  on  digitized 
tapes . 

The  stiffened  panel  components  of  an  aircraft  can  have  arbitrarily 
curved  cylindrical  geometry  and  nonideal  boundary  conditions.  To 
provide  better  approximate  techniques  for  analyzing  these  complex 
stiff etied  panel  structures,  N0VA-2S  was  modified  in  KA  TM-118*  and 
Reference  4  to  include  free  and  elastic  rotationally  constrained 


*Stagliano,  T.R.,  Mente,  L.J.  and  Lee,  W.N.,  "Memorandum  on  the  Support 
Work  for  the  STRESNO  Phase  II  Test  Program",  Kaman  AviDyne,  KA  TM-118, 
March  1979. 
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boundary  conditions,  discrete  linear  elastic  springs  and  forces,  and  a 
method  of  analysis  for  an  arbitrarily  curved  cylindrical  panel  by 
specifying  initial  radial  deviations  from  a  nominal  circular  shape. 
Thus,  any  combination  of  clamped,  pinned,  free  and  elastic  rotationally 
constrained  boundary  conditions  are  permitted  coupled  with  either  fixed 
or  free  inplane  conditions.  The  translational  springs  and  inplane 
forces  can  be  used  in  conjunction  with  the  normal  and  inplane  free 
boundary  conditions  to  represent  elastic  restrained  boundaries  with  or 
wi.hout  boundary  loading.  These  new  options  within  N0VA-2S  will  allow 
the  analyst  more  latitude  in  representing  the  actual  boundaries  of 
aircraft  panels. 


The  NOVA-2  response  results  have  been  compared  with  experimental 
data  from  unstiffened  and  stiffened  panels  that  have  well  defined 
geometry  and  boundary  conditions.  These  comparisons  are  presented  in 
Reference  5  and  show  very  good  correlation  between  analytical  and 
experimental  results  from  various  panels  loaded  by  pressure  pulses. 
These  results  gave  a  good  indication  that  the  theoretical  framework  of 
the  panel  analysis  in  NOVA-23  for  unstif fened  and  stiffened  panels  in 
both  the  linear  and  nonlinear  response  regimes  was  sound  when  the 
geometry  and  boundary  conditions  are  well  defined. 


The  objective  of  this  effort  is  to  determine  whether  the  K0VA-2S 
computer  code  is  applicable  for  aniyzing  actual  aircraft  stiffened 
panels  consisting  of  skin-stringer-frame  combinations.  Shock  tube  tests 
were  performed  oy  the  Boeing  Wichita  Company  in  1978  on  a  section  of  the 
KC-135A  aircraft's  fuselage  for  AFWL  and  on  the  aft  fuselage  and  bomb 
bay  doors  sections  of  the  B-52  aircraft  for  Defense  Nuclear  Agency 
(DNA) .  A  N0VA-2S  correlation  with  the  E-52  tests  results  is  given  in 
Reference  4  in  which  good  strain  response  comparisons  were  obtained  for 
the  B-52  aft  fuselage.  These  B-52  correlation  results  will  be 
contrasted  with  the  current  correlation  with  the  KC-135A  data  since  the 
common  two  types  of  fuselage  construction  me  -epresented;  namely,  the 
KC-i35A  fuselage  is  a  stringer-fuselage  svj  and  the  B-52  fuselage  is 
a  longeron-fuselage  system.  The  stringt.  system  of  the  KC-135A  fuselage 
is  characterized  by  numerous  stringers  attached  to  the  skin  and 
supported  by  floating  frames  (not  attached  directly  to  skin)  while  the 
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longeron  system  of  the  B-52  consists  of  four  longerons  supported  by 
frames  directly  attached  to  the  skin.  The  frame  spacing  in  the  K.C-135A 
fuselage  is  double  that  of  the  B-52  fuselage  which  is  a  normal 
difference  between  the  two  fuselage  systems.  Both  these  fuselage 
sections  provide  a  severe  test  of  N0VA-2S  capability  to  predict  the 
response  of  actual  large  aircraft  components  to  blast  overpressure 
loading. 

The  shock,  tube  test  program  for  the  KC-135A  fuselage  section  was 
used  to  evaluate  overpressure  damage  to  this  structure  with  and  without 
internal  pressurization.  Extensive  pressure  and  strain  experimental 
data  were  obtained  and  documented  in  Refer jnce  6  for  several 
overpressure  levels  at  two  blast  orientations.  The  measured  pressure 
data  from  these  tests  were  stored  on  digitized  tapes  which  are  used 
directly  in  this  effort  with  the  N0VA-2LTS  code  version. 

Section  II  of  this  report  presents  a  description  of  the  KC-I35A 
tests  and  the  overall  test  results.  The  N0VA-2LTS  analytical  models  for 
the  KC-I35A  fuselage  section  for  the  various  loading  cases  considered 
are  presented  in  Section  III.  In  addition,  new  modif  icatii.  .is  of  the 
N0VA-2S  computer  code  required  to  represent  the.  actual  structural 
behavior  are  presented  in  Section  III  and  the  Appendix.  The  comparison 
of  the  analytical  response  results  with  the  experimental  results  are 
given  in  Section  IV  for  the  seven  selected  cases  under  investigation. 
Section  V  presents  an  evaluation  of  the  N0VA-2S  (and  N0VA-2LTS)  computer 
code  for  predicting  the  structural  response  and  damage  occurring  in  the 
KC-135A  fuselage  section.  The  conclusions  and  recommendations  for  this 
effort  are  given  in  Section  Vi. 
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SECTION  II 

DESCRIPTION  OF  THE  KC-135A  TEST  PROGRAM  AND  RESULTS 


1.  GENERAL  TEST  PROCEDURES 

The  KC-135A  test  program  was  conducted  by  the  Boeing  Wichita 

Company  and  documented  in  Reference  6.  The  overpressure  blast  tests 

were  performed  at  the  THUNDERPIPE  Shock  Tube  at  Sandia  Corporation  in 

Albuquerque,  New  Mexico.  A  r-ft  section  of  the  KC-135A  fuselage  from 

BS  980  to  BS  1100  was  sealed  with  end  bulkheads  and  was  placed  in  the 

19-ft  shock  tube  and  subjected  to  generated  blast  waves.  The  upper  lobe 

compartment  of  the  fuselage  test  section  was  instrumented  with  pressure 

and  strain  gages  and  these  measurements  were  recorded  for  all  test 

conditions  covered.  Since  the  lower  lobe  of  the  fuselage  has 

intermediate  bulkheads  to  strengthen  this  compartment,  the  upper  lobe 

compartment  was  considered  the  structural  component  to  be  tested  and 

analyzed.  There  were  22  shots  in  the  test  series  at  the  Sandia  shock 

tube  at  two  blast  orientation  angles,  namely,  90  and  45°.  The  90° 

orientation  corresponds  to  the  blast  wave  intersecting  the  upper  lobe 

compartment  ut  the  center  of  the  crown  (overhead  burst  simulation)  while 

the  45°  orientation  simulates  a  burst  between  overhead  and  side-on. 

Eight  unpressurized  test  shots  were  performed  at  the  90°  orientation 

2 

from  overpressure  levels  from  0.3  to  3.5  lb/in  ,  and  two  unpressurized 
shots  were  performed  at  the  45°  orientation.  Twelve  pressurized  test 
shots  were  performed  at  the  two  orientations  at  various  combinations  of 
internal  pressure  and  incident  blast  overpressure.  The  upper  lobe 
compartment  of  the  KC-135A  test  section  was  pressurized  for  these  test 
shots.  From  these  22  shots  seven  were  selected  for  the  KC-135A 
correlation  with  N0VA-2I.TS  and  are  given  in  Table  1.  There  are  four 
test  events  at  the  90”  orientation  and  no  internal  pressure,  one  test 
event  at  90°  orientation  with  internal  pressure,  and  two  test  events  at 
45“  orientation  with  and  without  internal  pressure. 
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TABLE  1 


SELECTED  BLAST  SHL  S  FOR  CORRELATION 


Shot  No. 

Event  No. 

Blast 

Orientation 

(deg) 

Incident 

Overpressure 

(lb/in2) 

Internal 

Pressure 

(lb/in2) 

6 

78-112 

90 

1.52 

0.0 

17 

78-130 

90 

1.95 

0.0 

5 

78-105 

90 

2.50 

0.0 

22 

78-140 

90 

3.50 

0.0 

21 

78-139 

45 

2.10 

0.0 

16 

78-129 

90 

2.08 

2.5 

20 

78-138 

45 

1.85 

2.7 

2.  TEST  SPECIMEN  DESCRIPTION 

The  KC-135A  fuselage  is  a  stringer  system  type,  that  is,  a  large 
number  of  evenly  spaced  longitudinal  stringers  are  attached  to  the  skin 
and  supported  by  circumferential  frames  attached  to  the  inner  flanges  of 
the  stringers  as  shown  in  the  end  view  photograph  in  Figure  1  and  in  a 
more  detailed  view  in  Figure  2,  The  test  section  was  six  bays  long  from 
BS  860  to  B5  1]00.  End  bulkheads  were  attached  to  the  skin  ard 
stringers  of  the  fuselage  section  by  a  series  of  angle  clips  around  the 
circumference  of  the  section.  This  fuselage  section  consists  of  upper 
and  lower  lobe  compartments  separated  by  a  floor  system  between  WI.  202 


-r* 
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and  WL  210  as  illustrated  in  Figure  3.  The  upper  lobe  compartment  can 
be  pressurized  and  was  for  some  of  the  tests.  The  upper  lobe 
compartment  of  this  truncated  fuselage  section  is  the  component  that  was 
instrumented  to  measure  the  dynamic  response  from  a  simulated  blast 
overpressure.  It  should  be  noted  that  this  fuselage  section  was  part  of 
an  airplane  that  underwent  full-scale  fatigue  tests  prior  to  this  22 
shot  test  series.  Hence,  the  structure  of  the  test  specimen  had  been 
worked  and  certainly  represents  joint  conditions  common  to  those  in 
airplanes  that  have  logged  significant  flight  time. 

Figure  3  shows  two  views  of  the  KC-135A  fuselage  test  section. 
Figure  3b  illustrates  the  side  profile  view  of  the  upper  lobe 
compartment  while  a  cross  section  view  at  BS  1040  is  given  in  Figure  3a. 
In  the  six  bay  test  section,  five  Z-section  frames  of  7075-T6  aluminum 
are  spaced  20  in  apart  at  the  body  station  designated  in  Figure  3b. 

These  frames  are  attached  by  two  fasteners  (see  Fig.  2)  to  each  of  the 
3 1  hat  stringers  of  7075-T6  aluminum  which  are  evenly  spaced  around  the 
circumference  of  the  upper  lobe.  There  are  three  types  of  hat  stringers 
and  the  number  in  parenthesis  at  each  location  indicates  which  type 
exists  at  BS  1040.  It  should  be  noted  that,  in  general,  a  stringer  is 
not  homogeneous  over  the  entire  length  of  the  fuselage  section.  There 
are  splices  at  which  the  cross  section  changes  irom  one  type  to  another. 
The  cross  section  dimensions  of  the  Z  section  frame  and  the  hat  section 
stringer  are  given  in  Table  2.  The  type  of  skin  fastened  to  the 
stringers  varied  over  the  upper  lobe  Compartment  as  shown  in  Figure  3b, 
where  the  skin  is  divided  into  three  areas  indicated  by  the  dashed 
lines.  In  the  crown  area  above  stringer  S-7  the  skin  is  0.064-in 
7075-T6  aluminum.  Side  skin  of  2024-T3  aluminum  below  stringer  S-7  and 
between  BS  S80  and  1020  is  0.064-in  thick  while  between  BS  1020  and  1080 
the  skin  is  0. 051-in  thick.  Above  WL  230  the  fuselage  cross  section  is 
circular  with  a  radius  to  the  out- r  skin  surface  of  72  inches.  Below 
WL  230,  the  cross  section  deviates  slight j.y  from  a  circular  shape  and 
the  frames  are  deeper  in  this  region.  Above  WL  230  the  cross  sections 
of  the  frames  are  constant.  It  is  noted  that  the  frame  and  stringer 
cross  sections  are  formed  and,  therefore,  have  rounded  corners. 
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TABLE 


FRAME  AND  STRINGER  CROSS  SECTION  DIMENSIONS 


All  dimensions  are  in  inches 
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3. 


TEST  INSTRUMENTATION  AND  RESPONSE  RESULTS 


The  test  instrumentation  consists  of  pressure  transducers  and 
strain  gages  placed  on  the  upper  lobe  portion  of  the  KC-135A  fuselage 
test  section.  At  BS  1038,  22  pressure  transducers  were  mounted  to  the 
fuselage  at  approximately  10®  intervals.  Uniaxial  strain  gages  were 
placed  on  the  inner  and  outer  flanges  of  all  the  frames  of  the  test 
specimen  with  major  emphasis  on  the  central  frame  at  BS  1040.  Shear 
rosette  gages  were  also  mounted  on  the  web  of  the  frame  at  BS  1040.  The 
locations  of  all  this  instrumentation  are  given  in  Reference  6.  In  this 
correlation  effort,  emphasis  was  placed  on  the  response  of  the  flanges 
on  the  central  frame  at  BS  1040.  Thus,  Figure  4  shows  the  approximate 
locations  of  pressure  transducers  and  the  strain  gages  on  the  inner  and 
outer  flanges  of  the  central  frame.  The  inner  flange  gage  designated  by 
A  is  located  on  the  inboard  surface  and  the  outer  flange  ga^  designated 
by  B  is  located  on  the  outboard  surface.  The  rosette  gages  are  not 
shown  since  they  were  not  used  in  the  correlation.  The  more  exact 
positions  of  the  pressure  and  strain  gages  are  given  in  Table  3. 

The  pressure  data  for  each  shot  was  put  on  a  separate  magnetic  tape 
that  could  be  used  directly  by  the  NOVA-2  LTS  computer  code  Co  define 
the  loading  over  the  upper  lobe  fuselage  section.  These  loading  tapes 
were  used  to  create  data  files  which  were  accessed  through  the  AFWL 
computer  by  designating  the  event  number  of  the  shot.  The.  measured 
strain  time  histories  on  the  flanges  of  the  frame  at  BS  1040  are 
presented  in  plot  form  in  Reference  6.  The  pressure  time  histories  at 
each  position  are  also  given  in  Reference  6.  These  pressure  plots  were 
used  to  determine  if  any  gage  produced  bad  data,  so  thac  the 
corresponding  tape  channel  could  be  ignored  in  the  pressure  model.  It 
was  found  that  only  very  few  pressure  channels  had  to  be  eliminated  from 
all  the  data  sets  used  in  the  correlation  effort. 

4.  GENERAL  RESPONSE  BEHAVIOR  OF  FUSELAGE  TEST  SECTION 

The  test  shots  given  in  Table  1  were  selected  for  correlation 
because  they  represent  the  various  loading  conditions  covered  in  the 
test  series  presented  in  Reference  6.  In  the  test  series  the  test 
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O  PRESSURE  TRANSDUCERS  (BS  1038) 
—  STRAIN  GAGES  ON  FRAME  (BS  1040) 


TABLE  3 


PRESSURE  AND  STRAIN  GAGE  LOCATIONS 


Pressure  Gages 

Strain  Gages 

BS  1038 

BS  1040 

Gage  No. 

0-Position 

Gage  No. 

0-Position 

(deg) 

(deg) 

PS-1 

2.64 

G4-1 

3.8 

PS-2 

13.28 

G4-4 

21.0 

PS-3 

22.43 

G4-7 

37.1 

PS-4 

32.58 

G4-10 

50.9 

PS- 5 

42.33 

G4-13 

65.1 

PS-6 

53.07 

G4-16 

78.8 

PS-7 

63.02 

G4-19 

91-4 

PS-8 

73.16 

G4-22 

110.8 

PS-9 

82.91 

G4-25 

126.2 

rs-io 

92.02 

r>  /,  oO 

TOO  /. 
uu  •  *4 

PS-11 

102.86 

G4-31 

151.2 

PS-12 

113.0 

G4-34 

165.9 

PS-13 

123.91 

G4-37 

179.3 

PS-14 

132.95 

G4-40 

195.7 

PS-13 

142.6 

G4-43 

213.3 

PS-16 

152.85 

PS-17 

162.74 

PS-18 

173.48 

PS- 19 

182.89  . 

PS-20 

193.43 

PS-21 

202.68 

PS-22 

213.13 
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specimen  was  exposed  to  blast  overpressure  loadings  from  directly  above 
the  crown  and  45°  to  the  side  as  shown  in  Figure  3a.  The  longitudinal 
axis  of  the  test  specimen  was  always  parallel  to  the  blast  front.  For 
both  orientations,  four  level-  of  static  internal  pressure  (0.0,  2.6, 
5.9,  8.7)  were  investigated  to  determine  the  effect  of  internal 
pressurization.  Except  for  the  last  shot  (22)  ,  the  measured  strains  in 
the  frame  at  BS  1040  remained  in  the  elastic  range. 

For  the  blast  orientation  of  90°  and  no  internal  pressure,  incident 

2  2 

overpressure  levels  of  1.52  Ib/in  (shot  6),  1.95  lb/in  (shot  17) 

2 

and  2.5  lb/in  (shot  5)  produced  maximum  measured  compressive  strains 

of  3200,  5000  and  6600  uin/in,  respectively,  at  18.25°  to  the  right  of 

the  center  crown  line.  For  the  final  shot,  22,  at  3.5  lb/in^,  the 

maximum  inelastic  compressive  strain  is  estimated  between  40000  and 

50000  uin/in  at  16.55°  to  the  left  of  the  center  crown  line  before  the 

gage  failed.  Since  the  peak  response  is  significantly  off  center  in  all 

these  tests  and  has  switched  sides  during  the  final  test,  there  are 

probably  anomalies  in  both  the  symmetry  of  the  pressure  loading  and  the 

structure.  If  there  are  local  weak  points  in  the  structure  due  to 

fabrication  procedures  or  Che  prior  loading  history  of  the  fuselage, 

these  anomalies  cannot  be  analytically  modeled.  For  the  blast 

orientation  of  45v  to  the  right  of  the  center  crown  line  and  no  internal 

2 

pressure,  an  incident  overpressure  level  of  2.1  lb/in  (shot  21) 

produced  a  maximum  compression  strain  of  4200  uin/in  at  30.45°  to  the 

right  of  the  center  crown  line.  This  peak  response  position  is  14.55° 

to  the  left  of  the  loading  line.  It  should  be  noted  that  for  a  slightly 

2 

higher  incident  pressure  (2.1  versus  1.95  lb/in  ),  the  45°  orientation 
produced  a  lower  peak  strain  level  (-4200  versus  -  5000  uin/in)  than  the 
90°  orientation  loading.  When  internal  pressurization  was  used  in  the 
upper  lobe  compartment,  the  strain  levels  were  significantly  reduced  for 

the  same  incident  blast  overpressure.  For  example,  with  an  internal 

2  2 
pressure  of  2.5  lb/in  and  an  incident  overpressure  of  2.08  lb/in 

(shot  16),  the  peak  strains  ware  reduced  to  -1080  and  +1270  uin/in  from 

n 

the  -5000  uin/in  in  shot  17  (Ap  =  1.95  lb/in^)  with  no  internal 
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pressure.  The  strain  response  time  hictories  for  the  pressurization 

cases  exhibited  higher  frequency  oscillations.  Further  increases  in  the 

2  2 

internal  pressure  to  8.75  lb/in  from  2.5  Ib/in  only  reduces  the 

2 

peak,  strains  by  about  25%  for  the  2.0  lb/in  blast  overpressure  level. 

2 

In  the  final  shot,  22,  (Ap  -  3.5  lb/in  )  significant  permanent 
damage  occurred  to  the  frames  in  the  crown  region  of  the  upper  lobe 
fuselage  section.  A  detail  assessment  of  this  permanent  damage  is  given 
in  Reference  6.  Briefly,  in  the  crown  region  there  was  yielding  in  the 
outer  flange  of  the  central  frame.  At  several  local  positions  in  this 
crown  region  severe  yielding  of  the  frame  occurred  which  led  to  rupture 
of  the  frame  section.  This  severe  damage  was  precipitated  by  local 
plastic  buckling  of  the  outer  flange  of  the  frame  and  subsequent  plastic 
buckling  of  the  web  of  the  frame. 
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SECTION  III 

ANALYTICAL  N0VA-2LTS  MODELS  FOR  THE  KC- 135  FUSELAGE  SECTION 


Initially,  a  structural  model  of  the  KC-13SA  fuselage  test  section 
was  generated  assuming  that  the  pressure  loading  was  symmetric  for  the 
90°  orientation  case.  This  model  took  advantage  of  double  symmetry  and, 
therefore,  contained  16  longitudinal  hat  section  stiffeners  and  three 
circumferential  Z  section  frame  stiffeners.  This  first  trial  model  used 
34  modes  and  an  integration  net  of  19  by  17,  so  that  Central  Processing 
Unit  (CPU)  time  on  the  CRAY  computer  was  quite  reasonable.  From 
response  solutions  for  shot  5  u3ing  this  model  it  was  found  that  local 
buckling  of  the  skin  between  stringers  had  to  be  taken  into  account  in 
some  manner  that  would  not  require  a  prohibitive  number  of  modes  and 
integration  points  for  an  accurate  solution.  Furthermore,  it  was  found 
that  the  pressure  loading  for  the  90°  orientation  was  not  symmetric 
about  the  crown  line  and,  therefore,  a  full  model  in  the  circumferential 
direction  would  be  required  in  all  cases  for  an  accurate  solution. 
However,  the  symmetric  model  was  used  to  determine  respom  sensitivity 
to  various  parameters  such  as  spatial  loading  distribution,  boundary 
conditions,  buckling  criteria  and  skin  thickness.  Thus,  in  this  section 
the  general  structural  model,  the  skin  buckling  method,  the  response 
sensitivity  and  the  final  model  description  are  presented. 

1.  THE  NOVA-2I.TS  STRUCTURAL  MODEL  FOR  THE  KC-135A  FUSELAGE  TEST 
SECTION 

The  analytical  N0VA-2LTS  structural  model  for  the  upper  lobe 
compartment  of  the  KC-135A  fuselage  test  section  is  based  on  the  test 
specimen  description  given  in  Section  II--2.  For  the  analytical  model 
the  geometry,  boundary  conditions,  skin  thickness,  stiffener  cross 
sections  and  material  properties  must  be  specified.  For  the  actual 
structure  described  in  Section  11-2,  it  was  apparent  that  some 
compromises  would  be  necessary  in  the  modeling  for  N0VA-2LTS.  It  is 
recalled  that  the  crown  section  or  the  upper  lobe  compartment  is 
circular  above  WJ„  230  and  slightly  deviates  from  circular  between  WL  230 
and  210.  Since  the  major  strains  occurred  in  the  crown  region  of  the 


fuselage,  the  modeling  below  WL  230  could  be  compromised.  Hence,  the 
whole  upper  lobe  compartment  was  assumed  to  be  circular  with  a  radius  of 
72  in  to  the  outer  surface  of  the  skin.  The  subtended  angle  of  this 
stiffened  fuselage  panel  was  taken  as  215.9°  and  the  length  of  the  six 
bay  fuselage  test  section  is  120  in.  In  the  N0VA-2LTS  program  the  y 
coordinate  refers  to  the  1  ingitudinal  direction  and  the  6  coordinate 
refers  to  the  circumferential  direction.  In  the  y-direction  the 
boundaries  are  assumed  to  be  simply  supported  along  the  circumferential 
edges  that  are  attached  to  the  end  bulkhead  by  angle  clips.  In  the 
(S-direction  the  boundaries  are  assumed  to  be  clamped  where  the  frames 
intersect  the  floor  of  the  compartment.  Even  though  these  assumed 
boundary  conditions  are  probably  not  exact  and  could  be  modified  by 
allowing  elastic  torsional  constraint  at  the  edges,  the  boundaries  are 
far  enough  removed  from  the  crown  region  so  as  not  to  have  a  significant 
effect  on  the  major  strain  response ,  In  K0VA-2LTS  only  one  skin 
thickness  can  be  specified  over  the  entire  panel  while  from  Figure  2  it 
can  be  seen  that  the  thickness  varies  over  the  panel.  So  again  based  on 
the  fact  that  the  majui  strains  occurred  ir.  the  crown  region,  it  is 
assumed  that  the  skin  of  this  upper  lobe  compartment  panel  is  0.064-in 
7075-T6  aluminum.  The  reference  surface  of  the  panel  Is  set  in 
N0VA-2LTS  at  the  midsurface  of  the  skin,  so  that  the  radius  to  this 
reference  surface  Is  71.968  in. 

The  six  bay  test  specimen  panel  was  stiffened  by  five  internal 
frames  and  31  stringers.  The  frames  were  spaced  20  in  apart  and  the 
stringers  were  evenly  spaced  around  the  circumference.  In  the  final 
structural  model  symmetry  was  assumed  only  in  the  y  direction,  sc  that 
the  model  required  just  three  frames  to  be  specified.  The  hat  section 
stringers  are  attached  directly  to  the  skin  and  the  flanges  of  the  hat 
sections  are  attached  to  the  Z  section  frames.  Thus,  there  is  a  1.25-in 
gap  between  the  frames  and  the  lower  surface  of  the  skin.  The  stringers 
are  assumed  to  have  constant  cross  sections  over  the  length  of  panel 
based  upon  the  cross  sections  at  BS  1040.  The  type  of  cross  section 
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varies  around  the  circumference  as  indicated  in  Figure  2  and  Table  2. 
The  frames  are  all  the  same  for  the  structural  model.  Their  cross 
sections  are  constant  above  WL  230  and  variable  below  WL  230.  Both  the 
stringers  and  frames  are  formed  sections  which  contain  rounded  corners. 
In  defining  the  cross  sections  of  these  stiffeners  for  the  structural 
model  the  rounded  corners  are  approximately  taken  into  account. 

Table  4  gives  the  cross  section  input  dimensions  for  the  three  types  of 
hat  section  stringers  and  the  variable  7.  section  frame  used  in  the 
structural  mode.  For  the  variable  section  of  the  frames  the  cross 
sections  follow  those  given  in  Reference  6. 


The  stiffeners  and  the  skin  of  the  structural  model  are  7075-T6 
aluminum.  The  material  properties  used  in  the  model  are  as  follows: 


Modulus  of  Elasticity 
Shear  Modulus 
Poisson's  Ratio 
Mass  Density 

Yield  Stress  (elastic  model) 
Yield  Stress  (elastic-plastic 
Strain  Hardening  slope 


10.3  x  106  lb/in2 
4.0  x  106  lb/ in2 
0.3 

0.259  x  10~3  lb-s2/in4 
72000  lb /in2 
model)  76400  lb/in2 

1.85  x  105  lb/in2 


2.  SKIN  BUCKLING  CRITERION  METHOD 


From  ths  preliminary  model  of  the  KC-135  fuselage  section  for 
shot  5  it  was  found  that  if  the  skin  between  stringers  was  not  allowed 
to  buckle,  very  large  compression  stresses  were  generated  in  the  skin. 
These  0.064-in  skin  panels  between  stringers  are  separated  from  the 
frames  and  are  about  8.5  in  wide  in  the  circumferential  direction  and 
120  in  long.  They  will  buckle  elastically  in  the  circumferential 
direction  at  a  low  stress  level.  Therefore,  their  stiffness 
contribution  in  the  preliminary  model  was  much  too  high  and  this  caused 
a  very  low  response  level  of  the  fuselage  section  as  compared  to  the 
experimental  data.  Since  the  buckling  of  these  skin  panels  are  very 
important,  a  method  of  incorporating  this  effect  into  the  structural 
model  was  established. 
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TABLE  4 


STIFFENER  CROSS  SECTION  MODELS 


Stringers 


Segments  | 

Stringer 

1 

2 

3 

4 

Type 

h^(in) 

b^(in) 

h1(in) 

b^(in) 

h1(in) 

bi(in) 

h^in) 

bi(in) 

1 

0.051 

0.76 

1.0 

0.102 

1.199 

0.183 

1.25 

1.33 

2 

0.056 

0.935 

0.95 

0.112 

1.194 

0.201 

1.25 

1.58 

3 

0.064 

0.985 

0.95 

0.128 

1.186 

0.201 

1.25 

1.58 

Frames 


(h  *  1.25  in) 
o 


Segments 

1 

2 

3 

4 

5 

(deg) 

hi(in) 

b1(in) 

h .  (In) 
1 

bi(in) 

h^(in) 

b^in) 

( in) 

b^(in) 

( in) 

b^(in) 

0.0 

1.322 

0.872 

2.28 

0.144 

4.08 

0.072 

5.88 

0.072 

5.952 

0.87 

6.74687 

1.322 

0.872 

2.28 

0.144 

3.84 

0.072 

5.40 

0.072 

5.472 

0.87 

13.4937 

1.322 

0.872 

2.28 

0.144 

3.415 

0.072 

4.55 

0.072 

4.622 

0.87 

20.2406 

1.322 

0.875 

2.28 

0.144 

2.965 

0.072 

3.682 

0.075 

3.754 

0.87 

26.9875 

1.314 

0.8866 

1.58 

0.1104 

3.424 

0.064 

3.69 

0.1104 

3.754 

0.8866 

107.95 

1.314 

0.8866 

1.58 

0.1104 

3.424 

0.064 

3.69 

0.1103 

3.754 

0.8866 
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A  skin  buckling  criterion  method  has  been  developed  and  programmed 
into  N0VA-2S  and  N0VA-2LTS.  This  method  provides  a  practical  approximate 
solution  to  a  complicated  local  skin  buckling  response  that  occurs  in 
the  fuselage  section.  The  buckling  of  the  thin  skin  panels  between 
stringers  or/and  frames  of  a  fuselage  section  results  from  induced 
compressive  membrane  stresses  caused  by  the  engulfing  blast  pressure 
loading.  To  represent  the  local  skin  buckling  modal  patterns  directly 
into  a  N0VA-2S  structural  model  of  an  entire  fuselage  section  would 
require  a  prohibitive  number  of  integration  points  and  modes.  A  sn  ' ler 
time  increment  for  solution  would  also  be  required  which  would  result  in 
additional  computer  time.  It  should  be  noted  that  similar  problems 
would  also  exist  for  nonlinear  finite  element,  structural  models  in  that 
many  small  elements  would  be  required  to  represent  the  local  buckling. 

For  structures  in  which  the  response  of  the  frames  and  stringers 
are  of  prime  importance,  the  skin  buckling  criterion  method  offers  an 
attractive  alternative.  In  this  method  a  critical  buckling  stress  for 
the  skin  panels  is  selected  and  when  the  circumferential  compressive 
membrane  stress  in  the  panel  reaches  this  value,  the  membrane  stress  is 
reduced  in  some  prescribed  manner.  This  process  is  performed  at  each 
integration  point  on  the  skin  for  each  time  step  of  the  solution.  If 
the  stress  in  the  panel  becomes  lower  than  the  buckling  stress,  the  skin 
becomes  fully  active  again.  This  method  requires  the  definition  of  two 
parameters,  namely,  the  critical  buckling  stress  and  the  stress  decay 
constant.  The  selection  of  these  two  parameters  were  guided  by  static 
buckling  design  procedures.  Thus,  the  critical  buckling  stress  can  be 
computed  from  various  static  buckling  formulas  for  flat  and  curved 
panels.  It  should  be  noted  that  panel  buckling  is  influenced  by 
geometric  imperfections  (usually  unknown),  snap-through  action  (shallow 
curved  panels),  dynamic  effects  and  local  panel  boundary  conditions .The 
design  procedure  for  fuselage  sections  (see  Reference  7)  indicates  that 
the  effective  area  of  the  skin  panels  between  stringers  is  reduced  after 
buckling  by  the  ratio  of  the  critical  buckling  stress  (n  )  to  the 
applied  stress  (o  ).  Therefore,  the  stress  reduction  ratio  evaluated 
for  this  method  in  N0VA-2S  is  ( o  /o  )n  ,  where n  is  defined  as  the 
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the  stress  decay  constant.  For  the  8.5-in  x  120-in  panels  between 
stringers  of  the  KC-135A  fuselage  section,  the  static  critical  buckling 
stress  was  estimated  from  the  following  buckling  formula  for  a  clamped 
flat  panel  comp res sively  loaded  along  the  long  edges: 

it  k  E  ,h^2  ,,  ,,  2 

s  »  ■■  tt: (— )  —2100  lb  /in 

cr  12(,l~v2)  a 


where 

k  «  buckling  coefficient  *  4 

^  6  2 
E  »  Modulus  of  Elasticity  =  10.3  x  10  lb/in 

v  =«  Poisson's  Ratio  =  0.3 

h  =  skin  thickness  -  0.064  in 

a  <=  short  length  =  8.5  in 

In  the  final  NOVA-2LTS  models,  the  stress  decay  factor  (n)  was  set 
at  unity  for  the  best  correlation.  The  selection  of  n  °  1  is  consistent 
with  static  design  procedures  and  the  n  used  for  the  correlation  with 
the  B-52  aft  fuselage  section  data  in  Reference  4.  Thus,  n  was 
permanently  set  in  the  program  at  unity  and  it  was  not  introduced  as  an 
input  quantity. 

To  accomodate  the  skin  buckling  criterion  method  additions  in  the 
DEPROP  subroutine  of  the  NOVA-2S  and  N0VA-2LTS  versions,  the  following 
modifications  were  made  in  the  DEPROP  input. 

Group  10,  p.  1)6  of  Deference  3,  has  been  modified  as  follows: 

Group  10  (3112)  NL,  NSHEAR,  NBUCK 
Number  of  layers.  (NL) 

(NI.  must  be  1  for  KTYPE  =  1  or  2;  and  3  for  KTYPE  =  3  or  4  ) 

Core  shear  deformation  option  (KTYPE  =  3  or  4)  (NSHEAR) 

0,  no  shear  deformation 
1,  shear  deformation  included 
Skin  buckling  criterion  option  (NBUCK) 

0,  no  skin  buckling 
1,  skin  allowed  to  buckle 
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If  NSHEAR  ~  0  (Group  10)  skip  to  Group  17K,  unless  NBUCK  »  0 

The  following  data  group  is  included  in  Group  17. 

Skip  Group  17K  if  NBUCK  -  0 

Group  17K:  (F12.1)  BUCKS 

2 

Critical  Compressive  Buckling  Stress,  Ib/in  (BUCKS) 

(should  be  a  negative  number) 

The  changes  in  the  NOVA-2S  and  NOVA-2LTS  program  listings  to  extend 
DEPROP  to  include  the  aforementioned  addition  are  presented  in  the 
Appendix.  It  should  be  noted  that  the  core  shear  deformation  option  for 
honeycomb  panels  indicated  in  Group  10  is  an  in-house  addition  to  DEPROP 
that  Is  not  yet  complete.  Therefore,  users  should  .just  set  NSHEAR  to 
zero . 

Since  the  frames  of  the  KC-135A  fuselage  are  not  directly  attached 
to  the  skin,  the  skin  buckling  criterion  was  applied  to  all  integration 
points.  However,  in  the  case  of  the  B-52  aft  fuselage  in  Reference  4 
where  the  frames  are  directly  attached  to  the  skin,  integration  points 
along  the  frames  were  excluded  from  the  skin  buckling  option  in  oiuei  to 
allow  for  local  fully  effective  skin.  At  present  there  are  separate 
options  for  each  fuselage  model  which  have  been  handled  by  internal 
changes  in  the  program.  In  the  future,  the  fully  or  partially  effective 
skin  options  should  become  part  of  input  groupings  for  specified  liras 
of  integration  points. 

3.  STRUCTURAL  RESPONSE  SENSITIVITY  TO  VARIOUS  PARAMETERS 

In  the  preliminary  phase  of  this  effort  many  KC-135  fuselage  runs 
were  made  with  the  symmetry  model  at  relatively  low  CPU  time  to 
investigate  the  response  sensitit itv  to  various  parameters,  especially 
those  parameters  associated  with  the  skin  buckling  criterion  method.  A 
summary  of  all  these  runs  are  given  in  Table  5  which  was  organized  to 
indicate  the  response  sensitivity  to  the  spatial  loading  distribution, 
stress  decay  constant,  critical  buckling  stress,  boundary  conditions  and 
skin  thickness.  All  the  computer  results  shown  in  Table  5  are  for 
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TABLE  5 

KC-135A  FUSELAGE  STRUCTURAL  RESPONSE  SENSITIVITY 


(Shot  5,  Ap 


2.5  Ib/in2) 


Load  Model 

Boundary 

Conditions 

Code 

Stress  Decay 
Constant 

Buckling 
Stress  (psi) 

Skin 

Thickness 

(in) 

Maximum 
Strain 
(yin/ in) 

Loading  Distribution 

RS  Sym. 

21 

0.33 

-2000 

0.064 

■ 

LS  Sym. 

21 

0.33 

-2000 

0.064 

1 

Full 

21 

0.33 

-2000 

0.064 

jgm 

Stress  Decay  Constant 

RS  Sym. 

21 

0.0 

-2000 

0.064 

-5500 

RS  Sym. 

21 

0.2 

-2000 

0.064 

-6500 

RS  Sym. 

21 

0.33 

-2000 

0.064 

-6900 

RS  Sym. 

21 

0.5 

-2000 

0.064 

-7150 

RS  Sym. 

21 

1.0 

-2000 

0.064 

-7400 

i  Critical  Buckling  System 

I. 

Sym. 

21 

1.0 

-2000 

0.064 

-7400 

ym  ♦ 

21 

1.0 

-2500 

0.064 

-7500 

Rj  Sym. 

21 

1.0 

-3000 

0 . 064 

-7900 

RS  Sym . 

21 

1.0 

-3500 

0.064 

-7300 

rs  Sr  • 

21 

1.0 

-4000 

0.064 

-7800 

RS  S>  . 

21 

1.0 

00 

0.064 

-1500 

Boundary  Condition 

RS  Sym. 

11 

0.0 

-2000 

0.064 

J  1 

RS  Sym. 

21 

0.0 

-2000 

0.064 

PS  Sym. 

22 

0.0 

-2000 

0.064 

mm 

Skin  Thickness 

RS  Sym, 

21 

1.0 

-2000 

0.064 

RS  Sym. 

21 

1.0 

-2000 

0.0575 

ESI 
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shot  5  at  an  incident  overpressure  of  2.5  lb/ in  and  no  internal 
pressure.  The  load  model,  boundary  conditions  code,  stress  decay 
constant,  buckling  stress  and  skin  thickness  are  the  parameters  that 
were  varied  in  the  preliminary  analysis.  Table  5  shows  the  result  of 
varying  one  parameter  at  a  time  based  on  the  comparison  of  the  maximum 
strain  occurring  in  the  flanges  of  the  central  frame.  The  load  models 
considered  are  the  symmetry  structural  case  using  the  right  side  (RS) 
and  left  side  (LS)  loads  and  the  full  structural  and  load  case.  In  the 
boundary  condition  code  the  first  digit  indicates  the  condition  in  the 
y-direction  and  the  second  digit  the  condition  in  the  6-direction.  A 
number  one  (1)  designates  clamped-clamped  and  a  number  two  (2) 
designates  pinned-pinned. 

The  differences  in  the  spatial  loading  distributions  represented  by 
the  three  load  models  for  the  90“  orientation  in  Table  5  are  not 
significant  from  cursory  comparison  of  pressure  time  histories,  yet  the 
change  in  maximum  strain  response  can  be  as  much  as  27%  by  simplifying 
the  full  loading  case  by  assuming  symmetry.  If  the  stress  decay  con¬ 
stants  were  varied  from  the  linear  reduction  of  stress  (n  =  1)  case  to 
the  constant  stress  (n  =  0)  case,  the  maximum  strain  decreased  by  about 
26%.  It  was  also  found  that  If  the  skin  membrane  stress  was  allowed  to 
vanish  after  buckling  (n  *  «) ,  the  maximum  strain  increased  by  about  30% 

based  on  the  n  =  1  case.  When  the  critical  buckling  stress  was  varied 

2 

from  -2000  to  -4000  lb/in  ,  the  maximum  difference  in  strain  response 
produced  only  an  6.8%  change.  Thus,  the  maximum  strain  response  is  not 
very  sensitive  to  the  critical  buckling  stress  at  this  low  stress  level 
regime.  However,  when  skin  buckling  is  ignored  (o  =  00 ) ,  the  maximum 
strain  response  reduced  by  about  80%.  For  the  boundary  condition  sensi¬ 
tivity  evaluation  the  code  21  is  considered  the  base  case,  that  is,  pinned 
in  the  y-direction  and  clamped  in  the  6  direction.  If  the  boundaries 
are  clamped  in  both  direction,  the  maximum  strain  response  increased  by 
about  9%.  If  the  boundaries  are  pinned  in  both  directions,  the  maximum 
strain  response  decreased  by  about  4.5%.  Since  the  location  of  the 
maximum  strain  in  the  crown  region  is  far  removed  from  the  boundaries. 
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this  strain  response  was  not  very  sensitive  to  the  boundary  condition. 
Since  the  skin  thickness  varies  over  the  upper  fuselage  compartment, 
response  sensitivity  to  skin  thickness  was  considered.  A  structural 
model  using  an  average  skin  thickness  of  0.0575  in  increased  the 
maximum  strain  response  by  about  5.5%  over  the  model  based  on  the  crown 
region  skin  thickness  of  0.064  in. 

Based  partially  on  these  results,  the  final  N0VA-2LTS  models  for 

the  KC-135A  fuselage  test  section  used  the  full  loading  model  (ie,  no 

symmetry  about  center  of  crown),  a  stress  decay  constant  of  unity,  the 

2 

critical  buckling  stress  of  -2000  lb/in  ,  the  21  boundary  condition 
code  and  the  0.064-in  skin  thickness. 

4.  DESCRIPTION  OF  FINAL  STRUCTURAL  MODELS 

The  final  structural  N0VA-2LTS  models  for  the  KC-135A  fuselage 
test  section  were  constructed  based  on  the  information  generated  in 
Sections  III-l  through  III-3.  In  summary,  the  structural  model  uses 
symmetry  only  in  the  y  direction,  so  that  the  model  has  three  Z-section 
frame  stiffeners  in  the  ft  direction  and  31  hat  section  stringers 
stiffeners  in  the  y  direction.  A  uniform  0.064-in  skin  was  assumed  to 
e  attached  to  the  stringers  and  the  frames  are  attached  to  the  inner 
flanges  of  the  stringers  with  a  l. 25-in  gap  between  the  frames  and  the 
skin.  All  materials  are  assumed  to  be  7075-T6  aluminum.  The  fuselage 
test  section  is  120  in  long  in  the  y  direction  and  is  assumed  to  be 
circular  in  the  8  direction  with  a  radius  of  71.968  in  to  the 
midsurface  of  the  skin  and  a  subtended  angle  of  215.9°.  The  y 
boundaries  (circular  edges)  are  assumed  to  be  simply  supported  and  the 
8  boundaries  (straight  edges)  are  assumed  to  be  clamped.  The  inplane 
boundary  conditions  are  assumed  to  be  fixed  in  both  coordinate 
directions.  The  loading  model  used  the  appropriate  data  tape  and 
assumed  that  the  loading  distribution  in  the  y  direction  was  uniform. 

There  are  two  basic  N0VA-2LTS  models  that  are  repvesentive  of  the 
seven  shots  investigated  in  this  effort.  These  two  models  are  given  in 
Tables  6  and  7  and  are  based  on  an  elastic  response  model  (shot  5)  and 
an  elastic-plastic  model  (shot  22) ,  respectively.  The  data  deck 
listing  for  shot  5  given  in  Table  6  represents  the  model  used  for 
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TABLE  f> 

N0VA-2LTS  DATA  DECK  LISTING  THE  KC-L35A  FUSELAGE  F31  S-iQT 
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Table  6  (Continued) 
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Table  6  (Continued) 
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Table  6  (Continued) 
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elastic  solutions  of  the  90°  and  45°  orientation  loading  cases.  This 
basic  model  was  used  for  shots  6,  17,  5,  21,  16,  and  20,  except  for 
minor  changes  in  the  pressurization  cases.  This  model  employed  7 
symmetric  y  modes  and  20  6  modes  and  used  73  modal  combinations  for  the 
solution.  The  model  used  an  19  by  65  spatial  integration  net.  The 
time  increment  selected  for  solution  was  10  ms-  The  loading  input 
changes  for  the  various  shot  events.  The  start  times  for  shots  6,  17, 
5,  21,  16  and  20  are  0.3725,  0.3451,  0.3,  0.336/5,  0.3394,  and 
0.3458  s,  respectively.  In  shot  17  channel  13  (123.9°)  had  to  be 
eliminated  as  bad  data  and  in  shot  5  channel  3  (22.4°)  was  also 
eliminated  as  shown  in  Table  6.  All  the  pressure  data  for  the 
remaining  shots  were  fully  usable. 

For  the  unpressurized  shots  6,  17,  5,  and  21  the  critical  buckling 

2 

stress  was  set  at  -2000  lb/in  .  For  the  pressurized  shots  16  and  20, 
it  was  assumed  that  the  internal  pressure  would  prevent  the  skin  from 
buckling.  Therefore,  NBUCK  in  Group  10  was  set  to  zero  and  Group  17K 
(BUCKS)  was  eliminated.  In  order  to  include  the  internal  pressure  in 

shots  16  and  20,  KDS  was  set  at  3  in  Group  3  in  both  cases  and  Group  5 

2 

was  added  in  which  PS  was  sec  at  -2.5  and  -2.7  ib/'in  ,  respectively. 

The  data  deck  listing  for  shot  22  given  in  Table  7  represents  the 

KC-135A  fuselage  elastic-plastic  structural  model.  In  this  model  the 

modes,  integration  net  and  time  increment  are  the  same  as  the  previous 

model  for  shot  5.  Since  the  critical  buckling  stress  was  set  at 
2 

-2000  Ib/in  for  elastic  buckling  of  the  skin,  a  large  yield  stress 

value  was  inputted  to  keep  the  pseudo  skin  response  elastic  so  that  the 

overriding  buckling  criterion  would  not  be  confused  by  inelastic  skin 

response  quantities  in  the  program.  Therefore,  only  one  skin  layer  was 

necessary  in  the  model.  Since  the  stringers  remain  elastic  during  the 

response,  only  a  perfectly  plastic  material  model  was  inputted  with  the 

2 

yield  stress  at  72000  lb/ln  .  For  the  frames  which  went  inelastic,  a 

more  accurate  bilinear  representation  was  used  by  matching  an  actual 

stress-strain  curve  for  7075-T6  aluminum.  This  curve  fitting  resulted 

*> 

in  a  yield  stress  of  76400  lb/in*’  and  a  strain  hardening  slope  of 
5  2 

1.85  x  10  lb/ln  for  the  frame  material.  The  start  time  for 
shot  22  was  0.2949  s  and  all  pressure  channels  were  usable. 
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The  N0VA-2S  (N0VA-2LTS)  program  was  redimensioned  for  the  CRAY-1 
computer  at  AFWL  to  accomodate  the  models  given  in  Tables  6  and  7. 

These  dimension  changes  are  indicated  in  the  common  block  listing  given 
in  the  Appendix. 

The  preliminary  runs  for  shot  5  using  double  symmetiy,  42  modes 
and  a  19  x  33  integration  net  used  CPU  time  on  the  CRAY-1  computer  at  a 
rate  of  70  s/ms  of  response.  Some  runs  were  made  at  about  40  s/rcs  of 
response  by  reducing  the  integration  net  to  19  x  17.  The  final  runs 
using  73  modes,  a  19  by  65  integration  net  and  symmetry  in  only  one 
direction  used  CPU  time  at  a  rate  of  234  s/ms  of  response  for  shot  5 
and  292  s/ms  of  response  for  shot  22.  The  final  models  used  for 
correlation  are  probably  conservative  relative  to  general  accuracy.  By 
reducing  the  number  of  modes  to  about  60  and  0  integration  points  to 
33.  a  good  model  for  frame  response  could  be  run  at  a  rate  of  about 
100  CPU  s/ms  of  response.  It  should  be  noted  that  peak  response 
occurred  within  6  ms  of  response. 


SECTION  IV 

CORRELATION  OF  N0VA-2LTS  RESULTS  WITH  THE  TEST  DATA 


The  experimental  strain  results  of  the  central  frame  from  the  seven 
selected  test  shots  are  compared  with  the  corresponding  results  from 
NOVA-2LTS.  In  reviewing  all  the  data  it  became  apparent  that  there  were 
anomalies  in  the  actual  fuselage  structure  which  could  not  be  defined 
for  analytical  modeling.  Since  strain  response  is  very  sensitive  to 
irregularities  in  a  circular  curved  structure,  position  by  position 
comparisons  did  not  make  sense.  Hence,  for  each  test  shot  the  overall 
distribution  of  strain  maximums  around  the  circumference  of  the  upper 
lobe  compartment  of  the  KC-135A  fuselage  test  section  are  presented  from 
experimental  and  analytical  results.  In  addition,  the  strain  time 
histories  associated  with  the  largest  of  the  experimental  and  analytical 
strain  maximums  are  compared  for  the  outer  and  inner  flanges  of  the 
central  frame.  These  experimental  and  analytical  strain  time  histories 
do  not  necessarily  occur  at  the  same  circumferential  position.  The 
elastic  response  cases  at  the  90°  loading  orientation  (shots  6,  17,  5) 
are  presented  first  and  then  the  45°  loading  orientation  case  is 
presented.  Next  the  elastic-plastic  response  case  of  shot  22  is 
considered.  Finally,  the  correlation  results  of  the  pressurized  cases 
(shots  16  and  20)  are  presented. 

For  the  90°  loading  orientation  and  no  internal  pressure,  shots  6, 

17,  and  5  produced  elastic  response  of  the  central  frame  of  the  KC-135A 

upper  lobe  fuselage  section.  In  these  three  shots  the  fuselage  test 

section  was  subjected  to  engulfing  blast  waves  at  incident  overpressure 

2 

levels  of  1.52,  1.95,  and  2.5  lb/in  ,  respectively.  Figure  5  shows 
the  maximum  experimental  outer  flange  strain  of  the  center  frame  at 
three  key  positions  in  the  crown  region  of  the  fuselage  versus  the 
incident  overpressure.  These  three  measured  strains  from  gages  SG  19B, 
SG  22B,  and  SG  25B  indicate  that  the  largest  measured  strain  (SG  25) 
occurred  about  18  degrees  to  the  right  of  crown  center  line.  As  the 
incident  overpressure  increases  the  spread  among  these  three 
measurements  also  increases.  The  strain  value  for  the  measurement 
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Figure  5 


Maximum  Strain  Versus  Incident  Overpressure  for  90°  Loading 
Orientation  on  the  Unpressurized  Fuselage  Section 


SG  22B,  three  degrees  right  of  the  center  becomes  significantly  lover 

than  the  other  two  measurements.  This  drop  in  the  center  strain 

measurement  is  not  consistent  with  the  expected  strain  distribution  for 

a  circular  fuselage  section  or  the  analytically  determined  distribution. 

This  observation  was  the  first  indication  that  there  may  be  structural 

anomalies  in  the  crovn  region  of  the  fuselage.  Furthermore,  it  will  be 

shown  later  that  for  the  final  shot  22,  in  which  the  material  behavior 

was  inelastic  in  the  crown  region,  the  largest  strains  occurred  on  the 

left  side  of  the  crown  at  gage  positions  GS  13B,  GS  16B  and  GS  22B. 

Thus,  it  appears  that  the  experimental  strain  response  in  the  crown 

region  did  not  follow  a  consistent  deformation  pattern  as  the  level  of 

loading  was  increased.  The  maximum  strains  obtained  from  the  N0VA-2LTS 

analytical  model  are  also  plotted  in  Figure  5  for  incident  overpressures 

2 

of  1.52,  1.95,  and  2.5  lb/in  .  These  maximum  strains  occurred  at 

about  seven  degrees  to  the  right  of  the  crown  center  line  and  compare 

well  with  the  corresponding  maximum  experimental  values  of  gage  SG  25B 

(18  degrees  to  the  right  of  center).  Although  the  maximum  experimental 

and  analytical  strains  compare  very  well  there  are  differences  in  the 

strain  distribution  patterns  in  the  central  frame.  It  should  be  noted 

2 

that  shot  5  (Ap  *  2.5  lb/in  )  represents  a  deformation  condition  that 
is  near  the  threshold  of  yielding  for  the  frames  of  the  KC-135A 
fuselage. 

A  more  in  depth  understanding  of  the  strain  pattern  occurring 
experimentally  and  analytically  in  the  central  frame  is  obtained  by 
examining  the  results  generated  in  Figures  6  through  11  for  shots  6,  17, 
and  5.  In  Figures  6,  7,  and  8,  the  circumferential  distribution  of 
strain  maximums  for  both  the  outer  and  inner  flanges  of  the  central 
frame  are  given,  respectively,  for  shots  6,  17,  and  5.  In  these  plots, 
the  solid  and  dashed  Hues  represent  the  strain  distribution  from 
N0VA-2LTS  analytical  response  calculations  for  the  outer  and  inner 
flanges,  respectively.  The  circular  and  triangular  symbols  represent 
the  strain  distribution  from  the  various  strain  gage  measurements  for 
the  outer  and  inner  flanges,  respectively.  The  results  from  these 
comparisons  are  consistent  for  shots  6,  17,  and  5  and  the  following 
observations  are  made: 
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(1)  The  largest  compressive  strains  occur  on  the  outer  flange  and 
In  the  crown  region  of  the  frame  and  the  analytical  and 
experimental  values  compare  well  even  though  they  occur  at 
different  positions  about  12  to  16°  apart. 

(2)  In  the  crown  region  the  maximum  strains  distribution  predicted 
by  N0VA-2LTS  on  the  inner  flange  are  much  higher  (about 
double)  than  the  experimental  values.  However,  it  should  be 
noted  that  the  inner  flange  strains  would  not  be  a  governing 
factor  in  the  damage  criteria  used  in  N0VA-2S. 

(3)  The  irregularities  in  the  analytical  strain  distributions 
about  the  center  line  of  the  crown  are  caused  by  asymmetry  of 
the  loading,  while  the  experimental  distribution 
irregularities  may  also  be  attributed  in  part  to  structural, 
anomalies  as  well  as  the  loading  asymmetry. 

(4)  In  the  crown  region,  the  analytical  results  show  stronger 
membrane  action  across  the  frames  cross  section  than  the 
experimental  results. 

In  Figures  9  through  11,  experimental  and  analytical  strain  time 
histories  are  compared  for  the  outer  and  inner  flanges  of  the  central 
frame.  The  analytical  and  experimental  time  histories  are  given  for 
positions  at  which  the  largest  of  the  strain  maximums  for  each  shot 
occurred  in  the  previous  figures.  Experimentally,  the  peak  strain 
occurred  at  G  =  126.2°  (SG  25)  for  all  three  shots.  In  the  N0VA-2I,TS 
results,  the  peak  strain  occurred  at  S  «  114.7°  £tr  shots  5  and  17  and 
at  0  “  107.95°  for  shot  6.  The  general  shape  of  the  time  history  curves 
compare  well,  although  the  phasing  of  curves  differed  by  about  2  ms. 

The  maximum  compressive  strains  occurred  on  the  outer  flange  and  they 
compare  well  in  magnitude  for  all  three  elastic  cases  at  the  9C°  loading 
orientation.  Again  these  figures  demonstrate  that  the  N0VA-2LTS 
analytical  results  show  stronger  membrane  response  across  the  frame 
cross  section.  Possible  explanations  for  the  differences  in  the  strain 
distribution  across  the  cross  section  and  the  phasing  are  considered. 


67 


-4500 


-4000 


KC-I35A  FUSELAGE  CENTER  FRAME 

SHOT  6  £p  =  !.52L0/IN2 

- OUTER  FLANGEt 

r  N0VA-2LT 

—  • — INNER  FLANGE  J 
- OUTER  FLANGE  'l 


J-N0VA-2LTS  (0  =  107.95®) 


OUTER  FLANGE  'l 

INNER  FLAN6E}EXPER.MENTSG25(0  =  .26.2®) 


-3500 


-3000 


-2500 


-2000 


-1500 


-1000 


-500 


|\/)\ 
j  i  ' 


v  i 

ll  A 


TIME  (MS  ) 
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10.  Comparison  of  Peak  Strain  Time  Histories  at  the  90° 
Loading  Orientation  for  Shot  17 
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The  most  plausible  explanation  appears  to  be  out-of-roundness  of 
the  fuselage  section  in  che  crown  region.  Deviations  from  the  assumed 
circular  shape  would  produce  more  bending  action  around  the 
circumference  and  also  lower  the  frequency  of  response,  since  membrane 
response  is  of  higher  frequency  content  than  bending.  It  is  well  known 
that  the  responses  of  circular  shell  structures  are  sensitive  to 
out-of-roundness  (initial  radial  imperfections),  unfortunately,  the  true 
shape  of  the  fuselage  prior  to  each  test  shot  is  not  known.  If 
out-of-roundness  information  existed  for  this  KC-135A  fuselage  test 
section,  such  data  could  have  been  incorporated  into  the  N0VA-2LTS 
structural  model  through  che  initial  imperfection  option  in  the  code. 

The  inability  of  the  N0VA-2LTS  theoretical  formulation  to  properly 
predict  the  strain  distribution  across  the  frame's  cross  section  given 
the  correct  definition  of  the  actual  structure  was  discounted  because  of 
the  good  results  of  previous  correlation  efforts  in  References  4  and  5. 
In  particular,  the  correlation  of  the  B-52  aft  fuselage  results  showed 
that  the  frames  exhibited  very  strong  bending  action  in  the  region  of 
maximum  strain  which  was  predicted  very  well  by  N0VA-2LTS  when  the 
actual  geometry  of  the  B-52  fuselage  panel  section  was  modeled  through 
the  initial  imperfection  option. 

In  the  KC-135A  fuselage  the  frames  are  not  directly  attached  to  the 
skin  as  in  the  case  of  the  B-52  aft  fuselage.  The  N0VA-2LTS  formulation 
is  based  on  the  assumption  that  plane  sections  remain  plane  during 
deformation  of  the  whole  cross  section  which  includes  the  skin  and 
frames.  If  the  attachments  between  the  stringers  and  the  frames  became 
very  flexible  (sloppy)  during  the  response  this  assumption  could  be 
violated.  However,  since  the  skin  between  stringers  buckles  at  a  low 
stress  level,  the  effect  of  such  joint  flexibility  on  the  strain  and 
stress  distributions  across  the  floating  frames  should  not  be  an 
important  factor.  In  any  event,  such  undetermined  joint  flexibility 
would  be  difficult  to  model  accurately  with  any  structural  response 
code.  It  should  also  be  mentioned  that  errors  in  the  experimental  data 
are  not  uncommon,  in  particular,  phasing  and  polarity  errors  generated 


71 


during  the  data  acquisition  and/or  reduction  procedures.  When  the 
correlation  of  results  are  performed  soon  after  the  tests,  systematic 
differences  between  analytical  and  experimental  responses  can  be 
challenged  and  the  experimental  data  rechecked. 

Figures  12  and  13  illustrate  the  comparison  of  experimental  and 

analytical  results  for  shot  21  which  was  for  the  45°  loading  orientation 

2 

at  the  incident  overpressure  of  2.1  lb/in  .  Figure  12  shows  the 
distribution  of  strain  maximums  which  is  different  than  those  shown 
before  for  the  90°  loading  orientation.  The  region  of  the  largest 
strains  has  shifted  to  the  right  towards  the  direction  of  the  loading  at 
G  «  153°  for  both  the  experimental  and  analytical  results.  On  the  outer 
flange  of  the  center  frame  the  largest  of  the  strain  maximums  occurred 
at  Q  =  142°  and  155°  from  the  N0VA-2LTS  calculations,  while  the  largest 
experimental  strain  occurred  at  S  «  138°.  This  experimental  position 
favors  the  side  where  the  maximum  experimental  strain  occurred  for  the 
90°  loading  case  at  G  *  126°.  This  trend  again  indicates  the 
possibility  of  a  structural  anomaly  on  the  right  side  of  the  crown 
region.  The  largest  experimental  and  N0VA-2LTS  strains  differ  by  about 
10%  on  the  outer  flange.  On  the  inner  flange  in  the  region  of  maximum 
deformation  the  analytical  values  are  significantly  higher  than  the 
experimental  values  in  the  same  manner  as  the  90“  loading  cases, 
thereby,  indicating  more  bending  action  in  the  experimental  results. 
Overall  the  distribution  comparison  for  the  45°  loading  case  is  more 
consistant  than  the  90°  loading  case.  It  should  be  noted  that  with  the 
loading  direction  nearer  the  floor  support,  the  analytical  strains  at 
the  right  clamped  edge  are  higher.  Since  the  experimental  strains  at 
this  boundary  position  are  about  the  same  as  in  the  90°  loading  case 
(shot  17),  the  floor  joint  with  the  fuselage  is  probably  somewhat 
flexible.  Figure  13  compares  the  experimental  and  analytical  strain 
time  histories  of  the  largest  peak  strains.  The  phasing  differences  are 
less  for  the  45°  loading  case  than  the  90°  loading  case.  The  experi¬ 
mental  strain  distribution  across  the  frame's  cross  section  exhibits 
more  bending  action  than  the  analytical  results. 
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In  the  final  shot,  22,  (90®  orientation)  at  an  incident  over- 
2 

pressure  of  3.5  lb/ in  ,  the  frames  went  inelastic  over  a  large  region 
of  the  crown  and  at  four  positions  severe  plastic  buckling  of  the 
center  frame  occurred  leading  to  rupture  of  the  frame  at  three  of  the 
positions.  Figure  14  Illustrates  the  distribution  of  strain  maximums 
over  the  circumference  of  center  frame  for  the  outer  and  inner  flanges. 
Along  the  position  axis  cross  marks  indicate  the  locations  where  plastic 
buckling  of  the  frame  occurred.  The  experimental  strain  distribution  is 
misleading  in  that  the  severe  plastic  deformation  at  three  of  these 
positions  are  not  represented  on  the  plot  because  the  strain  gages  were 
outside  this  local  area  of  deformations.  At  9  *  91.4°  gage  GA  19  was 
located  in  one  of  the  local  buckling  areas  and  large  plastic  strains 
(over  4%)  were  measured  before  the  gage  became  inoperative.  The  three 
other  positions  (cross  marks)  would  have  shown  similar  plastic  response 
on  boch  flanges  if  gages  had  been  placed  there.  The  yielding  strain 
value  is  indicated  on  the  ordinate  axis  of  Figure  14.  The  experimental 
and  analytical  strain  distribution  are  different,  although  both  show  a 
12t-po  regi on  of  plastic  deformation,  from  45°  to  145°  experimentally  and 
from  60°  to  216°  analytically.  Experimentally,  the  zone  of  maximum 
response  is  shifted  more  to  the  left  side  of  the  crown  which  differs 
from  the  previous  elastic  shots  6,  17,  and  5.  The  experimental  results 
show  concentrated  locations  (plastic  buckling)  of  very  large  plastic 
strains,  (greater  than  4%)  while  the  analytical  results  show  a  more  even 
distribution  in  the  crown  region  reaching  about  2.5%  maximum  strain. 

The  N0VA-2LTS  model  does  not  at  present  allow  for  local  buckling  of  a 
stiffener,  hence  without  these  local  failures  present  in  the  response 
calculation  the  energy  is  absorbed  more  evenly  over  the  frames.  A 
method  to  account  for  local  elastic  or  plastic  buckling  of  stiffeners  in 
N0VA-2S  is  discussed  in  Section  V.  Fnr  there  strain  measurements  not 
affected  by  local  frame  buckling,  inelastic  strains  varied  between  0.8 
to  1.5%.  Analytically,  the  strain  varied  from  0.8  to  2.5%  in  the 
inelastic  region.  The  analytical  results  in  Figure  1^>  show  that  large 
strains  were  generated  at  the  right  s.de  clamped  position.  The 
experimental  values  near  this  position  are  still  quite  low,  as  they  were 
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in  the  elastic  cases.  This  indicates  that  for  this  large  deformation 
response  case,  the  clamped  condition  does  not  properly  model  this 
boundary.  It  appears  the  floor  system  does  not  provide  the  rigidity 
that  was  assumed  and  the  frames  may  be  primarily  acting  as  full  rings 
around  the  entire  fuselage.  To  remedy  this  situation  in  the  future  the 
entire  frame  could  be  modeled.  The  clamped  boundary  may  still  not 
significantly  affect  the  strain  response  in  the  crown  region.  Figure  15 
illustrates  selected  strain  time  histories  on  the  outer  flange  of  the 
center  frame.  In  this  figure  the  two  largest  measured  strains  are 
compared  with  the  two  largest  analytical  strains  in  the  crown  region. 

The  time  of  peak  plastic  strains  compares  better  than  in  the  elastic 
response  cases. 

In  shots  16  and  20  for  the  90°  and  45°  loading  orientations, 

respectively,  the  upper  lobe  compartment  of  the  fuselage  was  pressurized 

2 

at  2.5  and  2.7  lb/in  ,  respectively.  In  these  two  pressurized  shots 

2 

the  incident  overpressure  of  the  blast  was  at  about  the  2  lb/in  level 
which  is  comparable  to  the  unpressurized  shots  17  (90°  orientation)  and 
21  (45°  orientation).  The  pressurization  in  shots  1<S  *nd  20 
significantly  reduced  the  strain  response  in  the  frames  of  the  fuselage 
both  experimentally  and  analytically  as  predicted  through  NOVA-2LTS.  It 
is  recalled  that  for  the  unpressurized  cases  of  shot  17  and  21,  the 
maximum  measured  strains  were  at  the  -5000  and  -4200  pin/in  level. 

Table  8  gives  the  maximum  compression  and  tension  strains  measured  in 
shot  16  and  20  on  the  inner  and  outer  flange  of  the  center  frame.  It 
can  be  seen  that  these  strains  are  significantly  less  than  the 
corresponding  aforementioned  values.  Shot  16  gives  maximum  strain 
values  of  -1075  and  1100  pin/in  while  shot  20  gives  values  of  -1000  and 
1300  pin/in.  Table  8  also  compares  the  maximum  compression  and  tension 
strains  predicted  hy  N0VA-2LTS  with  the  corresponding  experimental 
values.  It  was  found  that  on  the  inner  flanges  for  both  shots  the 
N0VA-2LTS  analytical  strain  maximums  were  higher  than  the  experimental 
values,  while  on  the  outer  flanges  the  analytical  values  were  lower  than 
the  experimental  values.  Both  the  experimental  and  analytical  strain 
distributions  had  irregular  patterns  which  made  a  detailed  comparison  of 
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TABLE  8 


STRAIN  COMPARISONS  FOR  PRESSURIZED  KC-135A  FUSELAGE 
TEST  SECTION  IN  SHOTS  16  and  20 

Shot  16  -  90°  Loading  Orientation 
(ap  =  2.08  lb/in2,  *  2.5  lb/in2) 

i  Maximum  Frame  Strains  (yin/in) 


Position 

Experiment 

N0VA-2LTS 

Inner  Flange 

1100 

1740 

-900 

-1085 

Outer  Flange 

700 

540 

-1075 

-375 

Shot  20  -  45°  Loading  Orientation 
(Ap  *  1.85  lb/ in",  Pi  »  2.7  lb/in2) 


i  Maximum  Frame  Strains  (yin/in) 

Position 

Experiment 

NOVA-2LTS 

Inner  Flange 

1300 

1925 

-1000 

-1570 

Outer  Flange 

920 

395 

-800  I 

-464 

results  unfruitful.  The  importance  of  these  results  is  that  both  the 
experimental  and  analytical  strains  show  a  significant  reduction  in 
magnitude  from  the  internal  pressurization.  Many  of  the  differences  in 
the  results  may  be  the  result  of  unknown  structural  and  geometric 
anomalies  wnose  presence  were  suspected  from  previous  unpressurized  cor¬ 
relations.  Figures  16  and  17  show  comparisons  of  selected  analytical 
and  experimental  strain  time  hist-  les  for  shots  16  and  20, 
respectively.  In  Figure  16  for  the  90°  loading  orientation,  the 
analytical  strain  time  histories  at  the  two  positions  (0  =  87.7°  and 
67.5°)  of  maximum  compression  and  tension  strain  on  the  inner  flange  are 
compared  with  the  experimental  strain  time  history  at  the  position 
(0  -  110.8®}  for  which  both  the  maximum  compression  and  tension  strains 

occurred.  These  strain  responses  exhibit  higher  frequency  content  than 
found  in  the  unpressurized  response  cafes.  The  initial  strains  at  zero 
time  indicate  the  static  preblast  state  associated  with  only  the 
internal  pressure  loading.  In  Figure  17  for  the  45°  loading 
orientation,  the  analytical  strain  time  history  at  the  position 
(G  =  141.7°)  where  the  maximum  compression  strain  occurred  is  compared 
with  the  two  positions  (0  **  179.6°  and  195.7°)  at  which  the  maximum 
compression  and  tunsio-!  experimental  strains  occurred.  The  positions  of 
the  maximum  experimental  strain  are  well  away  from  blast  loading 
direction  of  0  =  153°,  while  the  analytical  strain  time  history  occurred 
near  the  loading  direction. 

Although  displacements  of  the  fuselage  test  section  were  not 
measured,  it  is  informative  to  present  some  of  the  radial  displacement 
results  from  the  N0VA-2LTS  response  calculations.  Figure  18  shows  the 
maximum  inward  radial  displacement  versus  ..he  incident  overpressure  for 
all  seven  test  case  considered  in  this  effort.  For  the  elastic 
unpressurized  case,  the  response  follows  a  linear  path  with  increasing 
overpressure,  but  becomes  nonlinear  as  inelastic  deformations  occur. 

The  displacement  for  the  elastic-plastic  response  case  of  shot  22  is 
substantially  higher  char,  that  of  shot  5  which  produced  response  near 
the  threshold  of  yielding.  The  maximum  displacements  of  the  pressurized 
cases  were  low  in  comparison  with  the  unpressurized  cases.  Figure  19 
illustrates  th?  differences  in  the  radial  displacement  distributions  at 
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the  time  of  maximum  displacement  between  the  90°  and  45e  loading 
orientations  for  shots  17  and  21.  Both  peak  displacements  are  very  near 
to  the  respective  positions  corresponding  to  blast  loading  direction. 
Figure  20  contrasts  the  radial  displacement  time  histories  for  an 
elastic  response  case  (shot  5)  and  the  elastic-plastic  response  case 
(shot  22).  Aside  from  the  peak  response  being  much  greater,  the 
elastic-plastic  response  case  has  a  longer  period  to  peak  than  the 
elastic  response  case. 
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SECTION  V 

EVALUATION  OF  THE  N0VA-2S  COMPUTER  CODE 

In  t.he  past  the  panel  subroutine  DEPROP  In  NOVA-2S  and  NOVA-2LTS 
has  been  evaluated  successfully  through  comparison  with  available 
experimental  results  and  other  nonlinear  computer  codes.  It  has  been 
established  that  NOVA-2S  can  accurately  determine  the  response  of  well 
defined  structures  within  the  framework  of  the  theoretical  formulation 
currently  programmed  in  DEPROP.  These  evaluations  have  been  well 
documented  in  References  2,  3,  4,  5,  and  8  for  unstiffened  and  stiffened 
panels.  Recently,  the  correlation  of  N0VA-2LTS  with  B-52  aft  fuselage 
results  from  blast  tests  performed  in  a  similar  manner  as  the  KC-135A 
fuselage  tests  has  been  completed  in  Reference  4.  These  results 
together  with  the  results  from  this  current  effort  are  now  employed  to 
assess  the  validity  of  using  the  N0VA-2S  stiffened  panel  option  to 
predict  the  response  of  large  complex  stiffened  aircraft  sections.  It 
should  be  noted  that  the  B-52  and  KC-135A  fuselage  sections  represent 
the  two  major  types  of  fuselage  construction,  namely,  the  longeron 
system  and  the  stringer  system,  respectively..  Thus,  the  correlation  of 
blast  test  data  from  these  two  types  of  fuselage  construction  with 
N0VA-2LTS  structural  model  response  results  covers  many  modelirg 
problems  chat  can  be  anticipated  for  actual  aircraft  sections.  It 
should  be  noted  that  some  of  these  modeling  problems  also  apply  to  any 
structural  code,  such  as  finite  element  code  representations. 

From  the  correlation  results  of  the  KC-135A  and  B-52  fuselage  test 
sections,  it  was  found  that  the  strain  response  of  the  frames  can  be 
significantly  influenced  by  the  boundary  conditions,  the  initial  panel 
geometry,  the  spatial  loading  distribution  and  the  local  buckling  of  the 
thin  skin  between  stringers  or  frames.  In  the  B-52  aft  fuselage 
correlation  in  Reference  4  the  N0VA-2LTS  structural  model  had  to  account 
for  the  noncircular  shape  of  the  fuselage  section  (through  the  initial 
geometric  imperfection  option),  local  skin  buckling  between  frames  and  a 
boundary  being  close  to  the  region  of  maximum  strain.  All  these 
ingredients  were  important  in  correlating  successfully  with  the 
experimental  data  using  the  modeling  techniques  provided  in  N0VA-2LTS. 
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In  this  current  KC-135A  fuselage  correlation  it  was  found  that  apparent 
minor  asymmetries  in  the  spatial  loading  distribution  for  the  9(J° 
orientation  were  important  and  therefore,  the  computer  time  advantage  of 
double  symmetry  could  not  be  used.  Furthermore,  local  skin  buckling 
between  stringers  was  very  important.  For  the  KC-135A  fuselage  the 
panel  geometry  and  the  boundary  conditions  were  not  major  factors  in 
influencing  reasonable  correlation.  The  designed  geometry  of  the 
fuselage  was  for  the  most  part  circular,  although  the  deformation 
henavior  pattern  indicated  that  perhaps  out-of-roundness  or  other 
structural  anomalies  were  present  in  this  fuselage  section.  Such 
possible  anomalies  could  not  be  defined  for  modeling  by  N0VA-2LTS,  but 
the  absence  of  finer  geometry  definition  did  not  prevent  a  reasonable 
correlation.  The  boundaries  of  the  KC-135A  fuselage  section  were  far 
enough  away  from  the  region  of  maximum  strains  that  the  assumed  boundary 
conditions  were  not  an  important  factor  for  the  elastic  response  cases. 
For  the  inelastic  response  case  the  assumed  clamped  condition  at  the 
floor  boundary  appeared  to  be  too  severe.  Whether  the  higher  local 
plastic  strains  generated  at  this  boundary  affected  the  strains 
calculated  in  the  crown  region  was  not  pursued  further,  since  the  more 
important  effect  of  local  plastic  buckling  of  the  frames  could  not 
presently  be  modeled  by  N0VA-2LTS.  N0VA-2LTS  should  be  modified  to 
account  for  local  elastic  or  plastic  buckling  of  stiffeners.  This  could 
be  accomplished  in  a  similar  manner  that  the  skin  buckling  was  handled 
in  N0VA-2LTS.  Each  cross  section  segment  of  a  stiffener  would  have  a 
critical  buckling  stress  parameter  modified  by  a  plasticity  factor  and  a 
stress  decay  factor,  sc  that  the  local  stiffness  of  each  segment  could 
be  altered  as  buckling  progresses  through  the  flanges  and  web  of  the 
stiffener.  This  is  a  recommended  modification  for  N0VA-2LTS  and 
NOVA-2 S. 

In  both  the  KC-135A  and  B-52  fuselage  testing  emphasis  was  placed 
on  measuring  the  response  of  the  frames  since  they  were  the  major 
structural  component  of  the  fuselage  panels.  The  thin  skins  of  the 
fuselages  were  undergoing  b,-ckllng  or  wrinkling  during  the  deformation 
of  the  stiffened  panel.  It  should  be  noted  that  once  the  skin  is 
allowed  to  buckle  in  N0VA-2I.TS  through  the  buckling  criterion  option  the 
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true  strain  response  of  the  buckled  skin  is  uncalculabie.  If  accurate 
strain  response  of  the  skin  was  desired  in  these  large  structurally 
stiffened  panels  a  very  fine  model,  probably  producing  a  prohibitive 
number  of  modes  and  integration  points,  would  result.  Whether  such  a 
solution  would  be  reasonable  on  the  computer  would  have  to  be 
determined.  The  same  problem  faces  the  finite  element  method  in  that  a 
prohibitive  number  of  nonlinear  finite  elements  are  probably  required  to 
simulate  the  local  buckling  patterns  between  stiffeners.  If  the 
computer  solution  is  unreasonable,  the  finite  element  method  would  have 
to  be  modified  by  special  birch  and  death  options  foi  the  elements  to 
approximately  account  for  skin  buckling  as  was  presently  done  in 
N0VA-2S.  It  should  be  noted  that  even  if  the  skin  responses  were 
directly  modeled,  the  undefined  initial  geometric  imperfections  of  the 
thin  curved  skin  would  probably  prevent  an  accurate  determination  of  the 
strain  response  of  the  buckled  skin. 

The  correlation  results  from  this  KC-175A  fuselage  investigation 
and  the  previous  B-52  aft  fuselage  effort  demoustrat  ?d  that  the  N0VA-2S 
or  N0VA-2LTS  computer  code  can  predict  the  peak  response  of  actual  large 
stiffened  aircraft  sections  with  reasonable  accuracy  (error  less  than 
10%)  in  the  elastic  response  region.  Thus,  for  the  threshold  of 
permanent  damage  criterion  in  N0VA-2S  based  on  the  onset  of  yielding, 
the  N0VA-2S  peak  strain  response  would  result  in  a  good  prediction,  in 
the  inelastic  response  region  leading  to  catastrophic  damage  based  on  a 
fracture  strain  criterion,  the  evaluation  is  more  complicated.  In  the 
B-52  aft  fuselage  tests,  the  frames  failed  under  tension  strain  on  the 
inner  flanges.  This  failure  appeared  (the  strain  gage  failure  at  about 
1%  strain)  to  be  premature  for  this  fuselage  test  section  and  resulted 
from  small  holes  that  were  present  in  the  inner  flanges  which  probably 
held  cable  harnesses  in  the  actual  aircraft.  One  longitudinal  set  of 
holes  in  the  inner  flanges  of  the  frames  was  very  near  the  position  of 
peak  tension  strain  -'nduced  by  the  blast  loading  and  all  fractures  on 
the  frames  occutre  '  .hrough  these  holes.  Thus,  the  fracturing  of  the 
frames  occurred  prematurely  as  the  result  of  stress  concentrations 
around  the  holes.  How  to  account  for  such  irregularities  in  the 
structure  (usually  unknown)  is  a  question  that  can  not  be  addressed  in 
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this  effort.  If  Che  holes  were  not  present  in  the  frame,  N0VA-2S  could 
give  a  good  prediction  for  the  B--52  fuselage  of  the  peak  inelastic 
response  into  the  strain  region  where  fracture  would  normally  occur 
through  a  tension  mechanism.  Experience  is  needed  in  selecting  the 
fracture  strain  for  the  catastrophic  criterion.  In  the  inelastic 
deformation  shot  for  the  KC-135A  fuselage  section,  the  severe  damage  to 
the  frames  was  caused  by  local  compressive  failures  (plastic  buckling 
of  the  frame)  which  lei  to  fractures  after  severe  buckling  deformations 
had  occurred  across  the  frame's  cross  section.  This  catastrophic 
compressive  failure  mechanism  cannot  presently  be  predicted  by  N0VA-2S. 
The  N0VA-2LTS  results  of  shot  22  showed  peak  compressive  strains  at  the 
2.5%  level  and  no  significant  tension  strains.  Therefore,  under  the 
present  catastrophic  criterion  based  or  tension  strains,  N0VA-2S  could 
not  predict  this  catastrophic  failure  of  the  KC-135A  fuselage.  Thus, 
there  is  a  real  need  for  a  compressive  catastrophic  failure  criterion  in 
N0VA-2S  along  with  a  method  to  permit  local  buckling  of  the  stiffeners 
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The  usual  clamped  and  simply  supporter,  boundary  conditions 
permitted  in  any  structural  response  code  can  only  approximate  the  true 
boundary  conditions  found  in  actual  large  aircraft  sections.  Farther 
refinement  can  be  achieved  by  introducing  o  her  flexibilities  at  the 
boundaries  through  translational  anu  torsional  springs  to  simulate  the 
stiffness  of  the  adjacent  s  ructure  outside  the  boundary.  In  this 
method  the  difficulty  is  es;  jmating  aprior.  the  proper  spring  constants. 
The  other  approach  would  entail  modeling  th  structure  beyond  the 
apparent  boundary  of  interest.  This  method  which  is  certainly  more 
accurate  has  a  significant  computer  time  penalty.  The  N0VA-2S  computer 
code  has  presently  the  capability  tc  pursue  all  of  these  options,  if 
required  to  obtain  reasonable  results.  Frcn  the  correlation  efforts 
involving  the  KC-135A  and  B-52  fuselages,  it  appears  that  good  results 
can  be  obtained  if  the  boundary  approa  hed  either  a  damped  or  pinned 
condition.  If  the  region  of  maximum  s  rain  is  far  removed  from  the 
boundary,  the  effect  of  the  boundary  is  minimized  for  elastic  response 
of  the  structure.  If  the  peak  response  is  occurring  near  a  boundary  it 
is  important  to  know  whether  the  actual  condition  approaches  clamped  or 
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pinned.  When  the  displacements  of  the  fuselage  sections  become  large 
when  inelastic  deformations  occur,  the  boundary  conditions  become  more 
significant,  but  more  evidence  is  needed  to  assess  the  overall 
importance  in  predicting  peak  response  away  from  the  boundary. 

In  analyzing  large  stiffened  panel  sections  from  aircraft 
uselajes,  it  is  important  to  be  able  to  handle  noncircular  geometry. 
With  the  use  of  the  initial  geometric  imperfection  option  in  N0VA-2S, 
arbitrary  geometry  ot  a  panel  can  be  approximated.  This  method  worked 
well  for  the  B-52  fuselage  panel  which  had  a  very  irregular  geometry  in 
the  circumferential  direction  as  shown  in  Reference  4.  The  method  for 
computing  the  necessary  radial  displacement  coefficients  for  a  surface 
which  deviates  from  a  specified  circular  shape  should  be  incorporated 
directly  in  the  N0VA-2S  code.  It  would  also  improve  the  geometry 
capability  of  N0VA-2S  to  introduce  a  conical  panel  option  that  would 
permit  better  modeling  of  a  nose  radome  type  of  structure  when  used  in 
conjunction  with  the  initial  geometric  imperfection  option.  The  conical 
shell  theoretical  formulation  already  exists  as  an  option  in  the  DEP1CS 
cylindrical  shell  code  from  which  DHPROP  was  originally  developed  es  a 
cylindrical  panel  code. 

With  the  increased  usage  of  composite  materials  in  aircraft  design, 
the  elastic  material  model  available  in  K0VA-2S  could  be  refined  to 
include  a  more  general  anisotropic  material.  The  present  orthotropic 
material  model  has  been  adequate  in  the  past  since  the  material 
properties  for  the  composite  material  are  usually  approximated  by 
specif  lying  the  material  in  terns  of  the  orthotropic  material  constants. 
A  more  accurate  description  of  a  composite  skin  would  contain  constants 
beyond  the  five  orthotropic  constants  and  they  could  be  computed  if  the 
detailed  laminar  construction  of  the  skin  is  specified.  It  would  also 
be  helpful  in  establishing  better  failure  criteria  for  composites  in 
N0VA-2S,  if  the  interlaminar  shear  stresses  would  also  be  computed  for 
the  composite. 

The  pressure  loadings  on  the  panels  generated  in  N0VA-2S  from  the 
nuclear  encounter  should  be  reviewed  and  modified  where  necessary  for 
the  varies  aircraft  components.  At  present,  for  panel  analysis  in 
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N0VA-2S  the  pressure  time  history  of  the  blast  loading  is  computed  at 
the  center  of  the  panel  and  the  distribution  is  assumed  to  be  uniform 
over  the  panel.  This  assumption  was  established  when  only  small  panel 
components  were  being  considered.  Now  for  large  stiffened  fuselage 
panel  sections,  the  actual  distribution  of  the  pressure  loading  over  the 
panel  becomes  important.  Thus,  a  new  spatial  loading  distribution  over 
the  large  panel  needs  to  be  defined  in  N0VA-2S. 


SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  general  conclusions  emanating  from  this  NOVA-2S  correlation 
effort  with  the  KC-135A  fuselage  test  results  are  delineated  as  follows: 

(1)  It  was  found  that  the  maximum  strain  response  on  the  center 
frame  in  the  crown  region  of  the  KC-135A  fuselage  was  most 
sensitive  to  the  local  skin  buckling  between  stringers  ard 
the  asymmetries  in  the  pressure  distribution  over  the 
fuselage's  surface.  Both  of  these  effects  were  taken  into 
account  in  the  N0VA-2L.TS  structural  and  loading  model.  The 
skin  buckling  effect  required  a  special  modification  to  the 
N0VA-2S  program  which  is  documented  in  this  report. 

(2)  In  the  new  skin  buckling  criterion  option  in  N0VA-2S  it  was 
found  that  the  peak  strain  response  was  more  sensitive  to 
the  stress  decay  rate  in  the  buckled  skin  than  the  critical 
buckling  stress. 

(3)  By  examining  the  experimental  strain  behavior,  it  is 
suspected  that  structural  anomalies  (most  likely  out-of- 
roundness)  were  present  in  the  crown  region  of  this 
fuselage  test  section  which  had  previously  beer,  subjected 
to  fatigue  testing.  Because  such  structural  anomalies  were 
not  defined  for  the  fuselage,  they  could  not  be  accounted 
for  in  the  N0VA-2LTS  model;  however,  the  assumed  perfect 
circular  geometric  shape  still  produced  reasonable 
analytical  and  experimental  strain  comparisons. 

(4)  The  main  difference  between  the  analytical  and  experimental 
results,  which  could  be  attributed  to  out-of-roundness ,  was 
that  the  experimental  strain  distribution  across  the  frames 
had  more  bending  content  than  the  analytical  results 
exhibited. 
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(5)  The  structural  model  employed  in  this  effort  represented  the 
upper  lobe  compartment  of  the  fuselage  where  the  boundaries 
at  the  floor  intersection  were  assumed  clamped.  This  model 
was  quite  adequate  for  the  elastic  response  cases,  especially 
since  the  maximum  strains  occurred  in  the  crown  region  far 
removed  from  these  boundaries.  However,  in  the  inelastic 
case  where  large  deformations  of  the  fuselage  occurred,  this 
clamped  condition  appeared  too  rigid.  It  appears  that  the 
floor  system  does  not  provide  full  rigidity,  since  the 

main  floor  beam  members  are  longitudinally  oriented.  The 
frames  are  continuous  through  this  floor  joint,  so  that  the 
assumed  clamped  condition  is  better  than  a  pinned  condition, 
but  there  is  an  undetermined  flexibility  at  this  boundary. 

This  flexibility  could  be  approximated  by  other  boundary 
options  in  NOVA-2S,  if  the  joint  stiffness  was  known  apriori. 
The  most  accurate  method  would  be  to  model  the  whole 
fuse  lage  and  absorb  the  computer  cost  penalty. 

(6)  N0VA-2LTS  predicted  the  maximum  strains  which  occurred  on 
the  outer  rlange  of  the  center  frame  very  well  (within  10£) 
for  the  unpressurized  elastic  response  case  at  the  90° 
loading  orientation.  For  the  45°  loading  case,  N0VA-2LTS 
predicted  the  change  in  strain  distribution  along  the 
circumference  and  the  maximum  strain  quite  well.  Fc-r  the 
pressurized  cases,  N0VA-2LTS  indicated  the  significant 
decrease  in  the  strain  level  that  has  occurred  experimentally. 
At  these  lower  strain  levels  the  comparisons  of  the  maximum 
did  not  show  the  accuracy  demonstrated  in  the  unpressurized 
cases.  This  is  not  unexpected  since  low  level  strain  response 
is  usually  more  sensitive  to  structural  anomalies. 

(7)  N0VA-2LTS  time  histories  of  the  maximum  strains  compared  well 
in  general  shape,  but  the  phasing  was  off  by  about  2  ms. 

The  experimental  time  to  peak  was  always  longer  which  could  be 
caused  by  the  greater  bending  content  across  the  cross  section 
of  che  frame  in  the  experimental  results.  Bending  is  a  lower 
frequency  response  than  membrane  action. 
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(8)  For  the  elastic  response  case  in  which  the  frames  failed, 
N0VA-2LTS  could  not  predict  the  local  plastic  buckling  of  the 
frames.  Analytically,  the  energy  was  absorbed  more  evenxy 
over  the  frames  and  the  compressive  strain  levels  reached 
about  2.5£,  while  experimentally  strain  levels  become  very 
high  at  local  positions  and  resulted  in  rupturing  of  the  frame 
at  several  circumfeiential  positions. 

(9)  Overall,  complex  aircraft  stiffened  panels  can  be  modeled  by 
N0VA-2S  or  N0VA-2LTS  with  confidence  provided  the  user  can 
define  all  important  structural  factors  involved  in  the  real 
structure  and  has  a  thorough  understanding  of  all  the  modeling 
options  now  available  in  N0VA-2S  and  how  to  apply  them  to 
represent  the  real  structure.  For  these  large  structures  the 
dimensions  of  the  program  had  to  be  extended  and  the  computer 
time  for  a  response  run  increased. 

Based  on  the  results  of  this  effort  the  following  recommendations 
are  given  as  follows: 

(1)  with  ali  the  modifications  introduced  into  NOVA-25  over  the 
last  few  years,  it  is  important  to  generate  a  new  User's 
Manual  that  contains  a  complete  group  listing,  gives  advice  on 
how  to  use  Che  new  modifications  and  gives  information  on  how 
to  redimension  the  program  for  modeling  the  larger  complex 
aircraft  components.  In  addition,  some  of  the  modifications 
require  that  changes  be  made  in  the  strain  and  stress  output 
format  of  the  program. 

(2)  NOVA-2S  should  be  modified  to  handle  local  buckling  of  the 
stiffeners  elastically  and  plastically.  In  addition,  the 
criteria  in  N0VA-2S  should  be  revised  for  catastrophic  damage 
to  include  a  compressive  type  of  failure  through  buckling. 

(3)  The  blast  induced  pressure  loading  formulation  used  in  N0VA-2S 
should  be  reviewed  and  modified  if  necessary,  especially  in 
specifying  a  non-uniform  pressure  distribution  for  large  panel 
structures.  There  are  other  less  important  modifications  that 
could  be  made  to  N0VA-2S  and  these  were  discussed  in 
Section  V. 
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(4)  Many  o£  the  new  additions  to  N0VA-2S  should  be  more  thoroughly 
checked  out  by  benchmark  solutions  using  other  structural 
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APPENDIX 


CHANGES  IN  I  HE  M0VA-2S  PROGRAM  L IS  I  IMG 


This  Appendix  gives  the  routines  in  NOVA-2S  and  NOVA-2LTS  that  were 
modified  for  increasing  the  program  dimensions  and  incorporating  the 
buckling  criterion  method.  The  COMMON  subroutine  in  the  Appendix  reflects 
the  dimension  changes  that  were  made  to  accommodate  the  structural  models 
generated  in  this  effort.  The  inclusion  of  the  buckling  criterion  method 
required  modifications  in  subroutines  DERV1,  DERV2,  DSET1,  DSET3,  SIGMA, 
and  SIGMAB  which  are  given  in  the  Appendix. 
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CnMM0<v/C8L  APK/  C010,Cim»CN8,09,EP80(1235)  .ETT.EXT.EXX, 

1  IA>Z(2)  ,KSUMA(  7095),  KY(  1235)  ,STT(  1235),  SXT(  1235)  ,SXX  (1235)  , 

2  Si  A (7095) ,S?a ( 7  0^5) ,S3A (7  095) ,S9 A (7095) ,S5A( 7095) ,S6A (7095) , 

3  uu(20, 13) , vv (20, 13) .WWC20, 1 3) . XK TT , XK X T , XKXX , X 1  A (7  095) , 

9  X2A  (7095)  ,X3A  (7095) .  X9  A  (  7095)  ,  X5A  (7  095) ,  X6A(7095)  , 

5  ZA(2).7b(2),ZF(b),ZG(6),ZH(fe) 

♦COkOECK  CBLK16 

COMMOM  /CP.UK  t  8/  A-EBGC31  ) .  AivE  HB  ( 1 29 ,  3 )  ,  C ,  CUR ,  F  T Z  1  (230)  ,FTZ2(230)  , 

1  FTZ3(230),FTZ9(230) ,FXZ1  ( 230 ) . FX Z2 (230 ) , FXZ3 ( 230 ) , F x Z9  (23u) , 

2  GB.GG.LrUT.LRHAR.LRG (10) ,L  R0UTC99) ,ISH (1 0) , LSH A R , LSOUT  ( 99 ) , 

3  LUSF.  (10,  10)  ,N3HEAR,PESTZ8(36J  ) , P£ ST ZG ( 3b 1 ) , PESX 78 ( 3b  1)  , 

9  PESXZG(3bl)  ,KSL(1 ') ,  1  0 )  , RSLP (  1  0 ,  ]  0)  .  RSM  ( 1  0 ,  1  0 )  ,  hS^P  (10,10), 

5  XLP7,XLP8»XJ>T2(3bl)»XSXZ(3b1  ) 

♦  DECK  0ERV1 

SUBROUTINE  0ERV1 

♦  CALL  C8LK1 

♦  CALL  CBLK2 
♦CALL  C8LK9 

♦  CALL  CHU1Q 

♦  CALL  CBLKU 
♦CALL  C6LM5 
♦CALL  CMOVA 
♦CALL  CBLANK 
♦CALL  CBLK16 

1  =  1 

on  100  P=1,A-G 

OU  10  0  Nsl,f'B 

T  F  ( MfSE  (iv  ,  R)  .  E  0 . 0  )  GO  TO  100 
IJU  (P,*)  =XX  (  I ) 
v  v  ( iv ,  6  ')  z  t  v  ( f;  f-:  B  + 1 ) 

*rr(N,M)=  /  X  f  1'G-' t:2+ 1  ) 

I  =  IM 
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101)  CONTINUE 
K  =  0 

00  500  Isl.MGT 
DO  500  J=1,NBT 

IFCNUSE(J.l).EO.O)  GO  TO  500 

K=Kfl 

SSlsO.O 

SS2=0.0 

SS3=0.0 

SS«=0.0 

SS5=0.0 

S$b  =  0.<) 

3S7=0.0 

SS8=C.O 

SS9=0.0 

SS10=0.0 

SSI  1  =  0.0 

SSI  2=0.0 

00  000  Msl.KS 

ms  (M-n*NGT*I 

Full  =  FUl(MM) 

FU22  =  FU2CNM) 

F V 1 1  =  F V 1 (NM) 

FV22  =  FV2(Niv.) 

T 1  =  FPl(NV) 

T2  =  FP?(NM) 

T3  =  FP3fr*) 

S1  =  0.0 

S3=0 . 0 

S0=0 . 0 

Sfe  =  0 . 0 

S7  =  0 . 0 

S9=0  . 0 

S 11 =0.0 

00  200  i*=t,NB 

IFCMUSF(N.M).E(J.O)  GO  TO  200 
NN  =  (N-l)  *fi§*T+J 
UMNsUU  r  N, V) 

VNnsVV  f  N , K ) 
fcPNs**  ( r -J ,  y) 

SI  =  I)MM*FU3(nM)  +  SI 
S3  =  UMn*FU«  Cmi\)  *  S3 
sa  =  vmn*fv3 (NNJ  ♦  su 

Sb  =  VPN*FVfl  (Ni'.j  +  Sb 
S7  =  ftN-f)*FP5  (Ni'j)  +  S7 
S9  =  *MN*FPb  (  I'm.)  +  S  9 
Sll  =  aiA"N*FP7  (ivtvj  +  Si  1 
200  CONTINUE 

551  =  S 1 *F  U 1  1  +  SSI 

552  =  SI  * F 0 ? 2  SS2 

553  =  S 3 * F l.1 1 1  +  SS  5 
SS«  =  S  u  *  F  V 1 1  +  SSti 
SS5  =  S  4  *  F  V  2  2  +  3  S  5 
SSb  =  S6*F  V 1 1  ♦  SSb 
SS7  =  S7 *  T 1  +  SS  7 
SSB  =  S 7  * T ?  +  oSB 
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SS9  =  S9*T1  +  SS9 
3310  s  S7*T3  ♦  SSI  0 
SSI  1  s  Sll *T1  +  SSl 1 
SSI?  s  S9 *  1 2  ♦  SS 1 2 
900  CONTINUE 
U(X)=SS1 
UG(K)=SS2 
UB(K)=SS3 
V(K)=SS4 
VG(K)=SS5 
V0(K)=SSf> 
w(K)=SS7 
WG(K)sSS8 
*6(K)=SS<> 

*GG(K)=SS1 0 
wBB  (K  J  =SS 1 1 
wG8lK)=SS12 
500  CONTINUE 

COMPUTE  STRAINS  AND  STRESSES 
Kso 

00  700  1=1. NGT 
K2  =  KSBX  f I ) 

00  700  .1  =  1  #  NBT 
IFCNUSE (J.I).EO.O)  GO  TO  70ft 
Kl  =  KSGXfJ) 

K  =  K  +  1 
UGF=UG(M 
UBFsUB  (« ) 

VF  =  V(M 
VGFsVG  (K I 
VBFsVB  ( K ) 

*F=W  (R  ) 
f.GFsiAG  (K) 

A8F S^B  (  K ) 

OwGF=0v\G  (R  ) 

OnfiF  =0^'6  (  R  ) 

EXX  =  X|_  1  *  f  UGF >XL1  *  (*GF*0WGF  +  0.5*(v'-GE**2-*’VGF**2  +  UGF**2)  )  ) 

ETTsXJ* (VBF+XJ* (AHF*OwflF+0.5*(*BF**2+VQF**2+UBF**?) ) ) 
EXT=XJ»N0F+XL1  *  IVbF*  ft  .0*XJ*VBF  )  ♦  X  J*  ( wGF  *  (rtHF  +  0  V.BF  )  +  D^GF  bF  +  UBF  * 
1UGF)  ) 

AC  =  *J*vRF  +  XLl*UGF  ♦  1.0 
*GGF  =  /;RG(*) 

W08F  = 

*G9F  =  *f,B  (*  ) 

XKXX  =  *L7  **GGF  *AC 
XKTT  =  X  J?*i-JftHF  **C 
XKXT  =  X!5**,jGhF*AC 
IF  (NPLT.FP.O)  bt-  I U  BOO 

E  T  T  =  E  T  r-.-.F  *  f  1  ,0  +  Xj*V6F-n#5*x-F)+VF*  (XJ*w8F  +  0.5*vF  T 
t!T=ETT  +  /F*e,*Uf  K) 

EXT  sEVT+XLI ★ ( AbF* VF-VGF*VF) 

XKTT  =  XKTT  +  XJ.J*vbF  +  XL1  *UGF-vjF 
XRXTsXKXT+XL 1  *  vbF 
AC  =  X  J  * -v E  F  +  -jF 
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XK XX  s  XKXX  -  XL'/*wGGF*wF 

XKTT  =  XKTT  ♦  X  J  *  VBF*  C  X  J*VBF  -  *F  )  -  XJ2*WF*WB8F  ♦ 
1  ( X J*  VBF  -  wF)<**2  ♦  AC*  (AC  *  *J**BF) 

XKXT  s  XK X T  -  XJb*WQBF*WF  +  XL3**GF*AC 
600  IF  (NSHEAR.EG.O)  GO  TO  605 
XI  =  XSX7 (K ) 

X2  s  XST7(K) 

EXX  s  EXX  +  X1*(0.5*X1  -  XL  1  * WGF ) 

£TT  s  ETT  +  X?*(0«5*X2  -  XJ*^BF  -  CUR*VF) 

EXT  s  EXT  -  XL1*"GF*X2  ♦  X1*(X2  -  XJ*^BF  -  CUR*VF) 

XKXX  =  XKXX  -  XL1*PESXZG(K) 

XKTT  =  XkTT  -  XJ*PEST28(K) 

XKXT  =  XKXT  -  XJ*PE3XZR(K)  -  XL  1 *PESTZG (K ) 

605  IF  (N0ERV.EG.2J  GU  TO  640 
CS1 1  =  CVJ 1 
OS  11  =  0*11 
FSl  1  =  FN’ll 
CS22  =  C^22 
DS22  =  D»'22 
FS22  =  FN22 


0S33=D«33 
SIG  =  1.0 

IF  (Kl.EQ.O)  GO  TO  610 
CS11  =  CM  1  +  C11GCK1) 

FS11  s  FMl  +  F11G(M) 

OSH  :  0^11  ♦  0 1 1 G  (K 11 
DS33=0S33  +  033G(*n 
SIG  =  STOEGIKIJ 
610  IF  (K2.FO.O)  GO  TO  620 
CS22sCMP?  +  C22e  (.J,*2) 

F322SFM22  +  F22B  ( J  ,t\2) 

0S2psr:«,2?  +  0226(J#K2) 

0S33=0 S 33+0338 (J ,K2) 

SIG  =  SIG+SIOEB (K2) 

620  S1A(K)  =  CS1  1  *Ex  X  ♦  CM 1 2*E  TT  +  FS11*XKXX  ♦  FM2*XKTr 
S2AIK]  =  CS22*£  T  T  +  CM2*FXX  +  FS22*XKTT  +  FMJ2*XKXX 
S3A(K)  =  CN'33*EXT  ♦  FM33  +  XKXT 
IF  (NBUCK  .  EG  .  0  .  OR  .  f-JC  ALL  .GT  .  0 )  GO  TO  625 
G2  =  CP22  +  ETT  +  CM2*EXX  +  Fm??*XKTT  +  FM2  +  XKXX 
IF  (G? . GT .BUCKS J  GO  TO  625 
S1A(K)  s  C  SI  1 *Ex  X  +  FS11*XKXX 

32 A (K)  =  f CS22-CP221*ETT+(FS22-FM2?) *XKT  T+BUCKS* (RUCKS /G2) **1 .00 
S3 A ( K )  =  CM33*£xT 
625  IF  (SIG. FG. 0.0)  GC  TO  630 

Sti  A  ( K  )  =  0S1  1  *XKXX  ♦  0M12  +  XKTT  ♦  FS!l*EXX  +  FM?*FTT 
S5A(K)  =  OS22*xkTT  ♦  0M2  +  XKXX  ♦  FS22*ETT  +  FM  2*EXX 

S6A(KJ=0S33*XkaT  +Fm33*FXT 

IF  (henCK.EQ.O.OK.NCALL.GT.O)  60  T 0  bb° 

IF  (G2.GT. BUCKS)  GO  TO  660 

SflA(K)  =  USH+xkxx  ♦  FSU+EXX  ♦  DM2  +  XKTT  +  F  M  2  *  E  T  T 
S5A(K)  =  0S22+XK I f  +  FS22+ETT  +  uM2*XKXX  ♦  FM1?*FXX 
SfcA(K)  =  0S33*XKXT  +  F  w  3  3  *  fc  X  T 
GU  TO  660 
530  SUA(K)  =  O.o 
S5A (K }  =  0.0 
S66(K)  =  0.0 
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GO  TO  660 


6«0  CALL  SIGMA  ( J  » I #  K  ) 

IF  (Kl.GT.O)  CALL  SIGMAB  (K1»0,J,I,K) 
IF  (K2.GT.0)  CALL  SIG^AB  fO»K2,J,I.K) 
660  X  1  A ( K J  =  EXX 
X2AOO  s  ETT 
X  3A (K )  =  EXT 
X  4 A (K )  =  XKXX 
X5AOO  s  XKTT 
X6A (K }  s  XKXT 
700  CONTINUE 
RETURN 
ENO 

♦  DECK  0ERV2 

SUBROUTINE  DFRV2 

♦CALL  CBLK1 
♦CALL  CRLK2 
♦CALL  CBLK3 
♦CALL  CPL*a 
♦CALL  CBLK10 
♦CALL  CBLK15 
♦CALL  C9LK17 
♦CALL  CNOVA 
♦CALL  C BLANK 
♦CALL  CBLK1P 


BOO  IZ=0 

DO  1B00  IPsl,RG 
MROsIIP-l )*NGT 

ou  lend  ts*i#«r 

IE (MuSECTS* 1 R) .EO.OJ  GO  TO  1600 

NNOs ( 1 S« 1 ) *  N  B  T 

IZ=I2+l 

SURSsO.C 

SVRSsO.o 

SHRSsO.O 

K  =  0 

00  1700  I=1,NGT 

MMSMMC+ I 

FUll  r  Fill  INN) 

FU22  =  FU2(N,v) 

F V  1  1  s  FV1  f 
FV22  =  FV?(Mm) 

T 1  r  FPlfl'M) 

T?  =  FF2  («■?•  ) 

T3  =  F  P  3  f  i  >  N ) 

SU  =  0. 

SV  =  0. 

S'*  s  n. 

pwlm  =  pr-’Acn 

on  thor.  .jsi.wRi 

IF  (NLSE  f  J,  I  )  .Fll-.o)  GO  TO  lfc<>" 

K  5  K  ♦  1 
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noo 


ir  tNuse r j. i) .eu.i )  so  to  itoo 

P«LN  a  PTNA(J) 

NNsNN(t+J 

UGFcUGOO 
U8F*UB (K) 

VFaV(K) 

VCFaVGOO 

VftFaVB(K) 

•iFiMfl' ) 

mGFcw6  (K) 

MdFs*B(K) 

WGGF  s  rtGGCK) 
wttBF  *  *BB(K) 

*GBF  s  «GB(K) 

0*GFa{)WG(K) 

0*8Fab*B(* ) 

!F(NU.FQ.O)  PPHxHm(K) 

PU  *  FU|l*Fu3CwU) 

PUG  *  Fll??*FU3lNNJ 
PUB  *  FU1  l*Ftia 
PV  a  FV1  l*FV3(U*v) 

PVG  a  FV??*FV3*NUJ 
PVB  a  FVn*FV«(Un) 

P*  a  Tl*FP5(NM) 

P*G  a  T?*FP5(imw) 

P«6  a  T 1  *FP6  ((*••«) 

P*GG  a  TT*FP5 (UU) 

PwBB  a  T 1  * F P 7 ( f . N j 
PWGB  a  T?*FPh('.‘wj 
PEXXUrXLl *PUG*C) ,U+XL1*UGF) 

PEXXVaxi7*VGF*HV6 
PEXXhsxi  7*P*G* («"GF ♦OwGF ) 

PFT  TUax  J2*l'BF*PuH 

PETTV  a  Xj*PVBMl.O  ♦  XJ*VBF) 

PETT*  s  Xj2*Pw8*  (*BF  ♦  HwfcF) 

PEXTUax  J*  (FUG*  U  .04-xt_l*l'GF)*XU  *uBF*PuGJ 
PFXTV  r  XL  1  *  (PtfG*  11,0  ♦  XJ*VBF)  ♦  X J*VGF*PVB ) 
PEXTIn  s  Xj3*(P*U*(«GF  +  r)fcGF)  ♦  P*G* C*BF+UrtBF ) ) 
PKXX*,s*L7  *PM>G 
PKTTWr  X  J?*Pv.H8 


PKXT^ax  J5*P<-GH 
AC  a  X J*VHF  ♦  XL1*L'GF 
PK  XXU  a  XL 1*XL /* nUGF*PUG 
PK XXV  a  X  J*XL  7*«;GGF*PVB 
PK XX1*!  a  PKKXw  ♦  XL  7 *P*GG* AC 
PkTTU  a  xJ2*XL1«v,bBF*PUG 
PKTTV  a  PVPoX  J*X 
PK T 1  *i  a  PkTT:<  ♦  P^BB^X J2* AC 
PKXTu  a  x ja*xL/»'*UBF*PUG 
PK  X  TV  s  XI  3*XJP*  •GHF*PVP 
PK  X  T  W  r  PKXTrt  t  X.JS*P*/GB#AC 
IF  (MSHF  Ak.E'i.o)  GO  TO  tlfcO 
XI  1  X  S  X  7  f  K  ) 

x?  a  xsrz(K) 

P E  X  X »■<  a  PF  <  X  -  -  xi.1*X1*PU(; 
PFTTV  =  PFTlV  -  CuN*X?*PV 
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■  PETTW  *  RE'TTft  -  XJ*X2*PfcB 
PpXTV  X T V  -  CUR*Xt*PV 

,f*EXT*  *^PEXTw  -  XL1*X2*PWG  -  XJ*XJ*PWB 

itad'sj  s  o# 

S2S0.Q: 

lF(Npt-T.EG.O)  GO  TO  1200 
PfcTTVsPFTTVtF V*i Vf+XJ*WBF) 

PETTtrisPETTt*  -  PW*  ( 1 .0  +  X  J*V&F-tvF)  ♦ 
PEXTVspf  XTV  ♦  XLl*(*GF*PV-^F*PVG) 
PEXTWsPFXTW  +  XL1*  CVF*PwG-VGF*Pft) 

PKTTU  s  XL1*PUG  +  PKTIU 
pkttv=xj«*pv&  ♦  pkttv 

PKTTWSPKTTW  -  P w 


X J**F*PVB 

X  J*  VF  *PWB 


PKXTV=XL1*PVG  ♦ 
PKXXft  s  PKXXft  - 
PKTTV  s  PKTTV  + 
1  PV*(3.*XJ*wBF 
PKTTlri  =  PKTT/i  ♦ 


PKXTV 

XL7* (PWGG*WF  ♦  Pft*ftGGF ) 
xJ*PV8*f4.*XJ*VBF  -  3.*wF)  + 

♦  2.*VF) 

X  J2*PX'8B*  f-WF  ♦  XL  1  *UGF )  ♦  PW*(2.*ftF  - 
1  3.*XJ*VRF  -  XJ2*w8BF)  ♦  P*R*X J*  f  a.*XJ**BF  +  3.*VF) 
PKXTV  =  PKXTV  ♦  XL3*PV*WGF 

PKXTw  s  PKXTft  -  x J5* (P^GB*HF  ♦  P  ft  *  <*  G  B  F  -  PftH**GF)  + 

1  PWG*XL3*  (XJ*i«HF  ♦  VF) 

51  s  Dwf)fK)  ♦  «F 

52  s  VF 

1200  PH  s  XLP?*PPP*PU* (aGF  ♦  DftGF) 

PV  s  XLP1  *PPP*Pv*  l  XJ*  (ftBF  +  HWirtF)  +  S2 ) 

Pft  =  XLP1 *PPP*Pn* (SI  -  XL  1 *UGF  -  X  J  *  vBF  -  1.0} 


FOX  s  o. 

F  V  X  =  0. 

FWX  s  0. 

SIGsl.O 

IF  (MnewV.LO.t )  GO  TO  1280 
IF  IKSTTF.E0.0J  GO  TU  1270 
ST IFFENFPS. 

K 1  =  KSGX(J1 
K2  =  ftSPXfl) 

IF  (Kl.FO.O)  GO  TiJ  1230 
SIGsSintGfM) 

IJ  S  (K1-1  1 

KSUM  S  KS»J*F?(IJ) 
wS  s  MSFGG (M  J 
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f  F  (K5UH.LT. "SJ  Gu  ro  1210 
ELASTIC  -  G  a 

G1  S  Cl  1  Gf  K  1  )  *X  1  A  (K  )  +  FH  G  (  K  11  *X<JA  (K  ) 

G«  =  f>UG(Kl)*A4AlK)  «•  F11G(K1)*X1  A  (K) 
G5=033G(«1)*aGA(KJ 

FOX=  (PEX  XU*51  ♦►•KX  x*J»GA  +  G5*Pk  XTM)  *C  A  1  +  FOX 
FVXS (PEX  X V*G l+PKx  X V*GU  +  G5*PK XTV)*CA1+FVX 
FwX=(PEXX'>*Gl  +  PKXX.i*G<J  +  G5*PKXT,v)*CA1+FWX 


GO  TO  1230 
tMF.LASTIC  -  G A P A  . 
1210  FSl  =  0. 

FS2  s  o. 

F  S3  *  o. 

1.0  =  fiSiwttX*(IJ-i) 
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00  1220  KKS1.N8 
L  =  L0  ♦  KK 

G!  =  AASGfKK,Kl)*SX(L) 

Z1  =  ZFRfKK.Kl) 

FSl  r  G1*(PEXXU  +  Z1*PKXXL)  ♦  FSl 

FS2  =  G1  *  (PEXXV  +  Z1*PKXXV)  ♦  FS2 

1220  FS3  =  Gl*fPEXXn  ♦  Z1*PKXXW)  +  FS3 

G2sCAHrl)33GfKU*XbA(K) 

FUX=FS1*CA2(K1 )+G2*PKXTU+FUX 
FVX=FS2*CA2(K1 J+G2*PKXTV+FVX 
FwXsFS3*CA2fKn+G2*PKXTW+FwX 

1230  IF  (K2.F0.0)  GO  TO  1270 
SIG=SIG*SIDEB(X2) 

IJ  s  MSGMP  +  (  K  2*  1 ) ♦NBAS  +  J 
KSIJM  =  KSUM0(IJJ 
NS  =  NSEGB(K2J 

IF  (ksijm. lt.ws)  go  to  1 2 a r» 

ELASTIC  -  BETA. 

G2=C229 ( J.K2) *X2A IK ) +F  228  t  J »  K2) *  X5 A ( K ) 
G5=D22B(J.*2)  *X5A  (K  )  ^228  ( J,K2)  *X2A  (K) 
G6=033B(J,K2) *AbA(K) 

FUX=(PETTli*G2  +  PK  T  rt.l*GS  +  PK  X  TU*Gb  )  *C  A  1  ♦K  UX 
FVX=(PETTV*G2  +  PKr  TV*G5+PKXTV*G8)*CA1+FVX 
FwXr(PETT'A*G2  +  PKTrw*G5tPKXTw*Gb)*CAl+FWX 
GO  TO  1270 
INELASTIC  -  BETA. 

12a0  FSl  s  0. 

FS2  =  0. 

FS3  =  0. 

KKKrwSG*  r K2- 1 ) +  J 
LO  =  f.'SMAX*(I.I-l  J 
OG  1 250  K  K: 1 , nS 
L  *  LD  ♦  kk 

G1sAAS9(kk,J,k2)*SX(L) 

ZlsZFO (KK.KKK) 

FSl  s  Gt*fPETTu  +  Zl*PKTTl.O  t-  FSl 

FS2  =  Gl*(H£TTv  ♦  Z1*PKTTV)  +  F S2 

1250  FS3  =  G1*(PETT.i  t  Z1*PKTTp)  +  F S 3 

G2  =  CAl*03  3H(Jf  K«f  J^XbA  (K) 
FGXsFSl*CA3(K?)i-G2*PKXTU  +  F(JX 
FVX  =  FS2*CA3(K2)»G2*PKXTV  +  F  VX 
FV' XsFS3*C  A3  ( *2 )  ♦  ts2*P*  X  TV*+F  vjX 


1270  J I  =  LPAP* (K-l j 
K  SUP  =  K S U M  A C  K  J 
IF  fKSliM.LT.LBAP )  GO  TO  1300 
IF  (KTfPF.EQ.3)  GO  TO  1300 


12*0  G 1 

2 

SI  A  f K) 

G2 

2 

S2A  f K ) 

G3 

z 

S  3  A  (  K  ) 

G« 

2 

SO  A  f K  J 

G5 

2 

ssa  r* i 

Gb 

2 

Sb A  f «  ) 

FI 

2 

»’F*Xi:*Gl 

♦  Kt  mi*G2  * 

PE  X ru*G3 

F  2 

2 

PK<<n*Ga 

+  pk  rro*G5  ♦ 

PK  X  r»J*G6 
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F3  s  PEXXV*G1  ♦  PETTV*G2  ♦  PEXTV*G3 

Fa  s  P*XXV*Ga  ♦  P*TTV*G5  ♦  PKXTV*Gfe 

F 5  s  PEXXW*G1  ♦  PtT  Tl^*G2  ♦  PEXTh*G3 

F6  s  PKXX**Ga  ♦  P*rTw*G5  +  PKXT**G6 

FU  s  CM  0*F  1  ♦  CM1*F2 
FV  s  CM0*F3  ♦  OU*Fa 
FW  s  CN 1 0 *F5  ♦  C  n  1 1  *  F  6 
GO  TO  1500 

1300  TOTUMso.O 
TOTVMs<t.f. 

totub  =  o.o  Reproduced  From 

TOTvb  =  o.o  Best  Available  Copy 

T0T*8  =  o.o 
DO  laoo  KK=l,LbAk 
L  -  J I  +  K  K 
SI  =  HGOfKK) 

S?  =  GX (KK ) *S1 

G 1  =  S X X  f(.  ) 

G2  s  STT(D 

G3  =  Sxr(L) 

TOTUM  =  TOTUr*  +  Sl*(PEXXU*Gl  ♦  PETTU*G2  +  PEXTu*G3) 

TOrUB  =  TOT  HH  t  S2*  t  PK  X  XiJ*G  1  ♦  PMTn*G2  +  PKX7U*G3) 

TOTVM  =  T(iTV»  ♦  5l*(PEXXV*Gl  ♦  PfcTTV*G2  +  PExrv*G3) 

TOTVB  =  TOTV8  ♦  S2*IPKXXV*G1  ♦  PKTTV*G2  +  PKXTV*G3) 

TOTWM  =  TOTrtM  ♦  Sl*(PEXX**Gl  ♦  PETTa*G2  ♦  PEXT.v*G3) 

MOO  TOTwB=TOT‘«e  +  S2*  (PKXXrt*Gl  +  PKTT,**G2  *  PK  X  T  v*G3 ) 

FU  =  CN**T0TUK  ♦  CNS* TOTUB*SIG 

F v  =  ro^*TUTVM  ♦  C‘19*TnTVn#.SIG 

F*«  =  CN«*  TOTr.M  +  CN9* TOTrtH *S1  G 


l  50C 

SU  = 

(Fll  ♦  P'J 

♦  FuXl*PRLN 

♦ 

SU 

SV  = 

TFV  ♦  FV 

♦  FVx)*PWL'N- 

* 

•SV 

S'*  = 

( F  *i  ♦  H 

*  F  >  X  i  *PFL'V 

■* 

SW 

IbOO 

CONTI 

HE 

SUkS 

=  SiikS  ♦ 

P  K  L  • -1  *  S  u 

SVKS 

=  S  V  R  S  ♦ 

PPL  ••1*5  v 

S*PS 

=  SvPS  ♦ 

1700  CONTINUE 

ELASTIC  SPPINGo. 

IF  (VSPP.EO.n)  go  TO  170" 

DO  l  730  II?  :  1,-vSP* 

LI  =  .\SPGfLl?) 

L2  =  NSPML12) 

K  =  (Ll-I  1»'HT  +  t_2 
BTGKL  =  *IGK  (1.12)  *<L?3 
IF  CIuIk(L12J  -  2)  1710, 1720. 17?5 
i  7 1 0  SUkS  =  Rll,KL*i»(i»)*Fui(MMO*i.n*Fu3(NNO  +  L2)  +  SUPS 
GO  TO  1730 

1  7  ?0  SVPS  =  a  IGAL*V  (t\  J  *FVl  (MMO-t-L  1  )  *Fv3  f  MNO  +  L2J  ♦  S  VPS 
GO  TO  1730 

1/25  S^'PS  =  RTf:kl*MM  *FPi  fwMOH.  1  >  *F  P5  ffj'M.H-L2)  *  S -PS 
17  30  C 0 rj T  I 1 1 E 

T  i )  W  S I ! !  ?  j  A  L  SP*  I-'-GS. 

17  40  IF  ( ->T  SC  .F ,J . •! )  :><J  TO  17  50 
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on  1746  i  =  i, ntsc  Reproduced  From 

J  =  TSLOCfl)  Best  Available  Copy 

II  =  KTSGB { I ) 

GO  TO  (1742,1743,1744,1745),  J 

1742  K  =  II 

SwRS  =  SW»S+XLP4*9IGC(I)**vG(K1  ftFP2(MM0+l)*FP5(NN0+II)*HK  (II) 

GO  TO  1746 

1743  K  =  (NGT-1 )*N6T  ♦  II 

SWRS  =  SwRS+XLP4»SIGC(I)*tNG(K)  *FP2(MMU+NGT)  *FP5(N('I0+II )  *HK(II) 
GO  TO  1746 

1744  K  s  (II-1)*NBT  ♦  1 

SWRS  =  SWRS+XLP5*8IGC(T)*wB(K)*FPl  (f^O+I  I )  *FP6  (NNO+l  )  *HJ  (  1 1 ) 

GO  TO  1746 

1745  K  =  II*N8T 

S*RS  =  SnRS*XLP5*8IGC(I)*nBCK)*FP1  (mn«o*I  I )  *FPb  ( MM0  +  N«T  1  *hj  (I  I ) 

1746  CONTIUflE 

1750  IF  (NIPF.EQ.O)  GO  TO  1780 
00  1770  1=1,4 

IF  (NFI(I).EQ.O)  30  TO  l?7o 
IN-PLAMF  EXTERNAL  FORCES. 

CALL  IHFORC  (F IP ( 1 . I ) . TIP ( 1 , I) , JL IP ( I) , NF I ( I) , T 1 *E , F ) 

IF  (1-2)  1752,1754,1760 
1752  F  =  F«XLP6*FU1 («hOM  ) 

GO  TO  1756 

1754  F  =  -F*XLPb*FUl MHO+NGT) 

1756  00  1758  J= 1 , NH  f 

1758  SURS  =  SOWS  -  F *SF IP ( J . I ) *FU3 (600+ J) 

GO  TO  1770 

1760  IF  (I. £0.41  GO  Tu  1762 
F  =  F»XLP6*FV3(  <rjU+l  ) 

GO  TO  1764 

1762  F  =  -F*XLP6*Fv3(i.NU+’JBT) 

1764  00  1766  .1  =  1,  NG  T 

1766  SVKS  =  SVRS  -  F*SFIP(J,T)*FV1 (MMO+J) 

1770  CONTroUfc 

1780  IF  (AOSCSsPS) .3T.1.0E20)  GO  TO  2150 
IF  (KCOUP.GT.o)  GO  TO  1790 
rr(IZ)  =  S'.'RS*PRL'JIIZ) 

Yf(TZ4.MGMR)  =  oVKS*PRLV(IZ) 

Y  Y  ( IZ+f-:Gfc,H2)  =  a.«n5*KKLWC  I Z  > 

GO  TO  t*0o 

1790  FHS(IZ)  =  SU6S 

HHS(  I7.+^G‘*6  J  =  Sl/i<  s 

kmS(  IZ+“G”02)  = 

1800  CONTINUE 

IF  (KCOUP.F.a.n;  GO  TO  2200 
DO  1 Ro 0  Mri.MG^H 
01  =  0. 

02  =  0. 

03  =  <•„ 

no  1850  G=i,f*G-«U 

01  =  K  '4  s  ( i-i )  *  i,  ij  i  n  ,  v.  j  +  r>  1 

02  =  RHS(fJ  +  VGtJ-j*AV(.>l,»*)  +  D? 

1850  03  =  C  H  S  f  ”  ♦ '  ’G  ^  H  2 )  *  A  *  ( N  ,  M  )  +  03 

Y  Y  CM)  =  i)t 
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Y Y (M+MGWR )  =  02 

1900  yy  (MfMRMBa)  =  u3  Reproduced  From 

go  TO  2200  Best  Available  Copy 

2150  KERR  =  1 
WRITE  (b, ? 151 } 

2151  FORMAT  (30HOS0LUUON  DIVERGING  IN  DEPROP  ) 

2200  RETURN 
END 

♦DECK  0SET1 

SUBROUTINE  USETt 

♦  CALL  CBLK1 

♦  CALL  CBLK2 
•CALL  CBLK3 
♦CALL  CRLKO 
♦CALL  CBLK5 
•CALL  CHLK7 
♦CALL  CBLK9 

♦  CALL  CBLMO 

♦  CALL  CBLKU 
♦CALL  C6LK13 
•CALL  CBLK15 
•CaLL  CNOV1 
♦CALL  CBLANK 
•CALL  CBLK10 

INPUT  DATA 

READ  (S.  70^01  MG#i,i1#MHAP,NbAP,LBAR 
READ  (5.7000)  (UG’UI)  .  1  =  1, f-G) 

REAP  (5.7000)  (wUNCI ) . 1=1 ,*6) 

READ  (5.7000)  NSY:<G,mSYM8,nHUV 
IF  (nKUV.EN.O)  t.b.iv  =  11 

READ  (5,7000)  PPL  r , NHnD . NPFP V , NSC . NSH , KCOUP 
READ  (5, 700  0)  r  „('ul 
IF (NNDUT.ED.O)  Go  TO  70 
UC  SO  1  =  1,  I 

50  READ  (5. 7000)  fUJUT  ( I ) ,  nOUT  1 1) 

70  READ  (5. 7(  oo)  i^P 

IF  (MkP.FD.o)  Go  TU  90 
DO  8  0  1  =  1  ,r-KP 

80  wfcA0(5,7O0O)  KPG ( l ) , KP8 f  I  ) 

90  IF  (KDAf'.FD.  !  .  AOO.MYPt'.EQ.  1  )  i-'1.'EkV  =  2 

IF  (KOAP.EU.  1  ..aoD.KTYPE.EN.3)  M'EPV  =  2 
IK  (M EPV.F','. 1 )  LEAR  =  I 
READ  (5,7000)  \'L,  I'-Sri&AR.NBl  C* 

IF  C (K TYPF+1 j /2.; E.2)  NSHEAR  =  0 
LR6AR  =  1 

IF  (KTYPt.LT.5j  f.L  =  3 
IF  (KTYpF.LT.3i  :<l  =  1 
RE  AO (5, 7 1  on)  XLP.IhETAO.A 
IF (NPL1 .FQ.o)  m=1.0 
IF  (Ni'E»-'''.EN.1  )  GU  IU  150 
f»P  1  Or  i  =  i  ,i  L 
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luo  REAP(5,7100)  HM(I),RHOM(I),EM(I)  ReDr< 

READ  (5#  7100)  TM'.SIGO.EP.EPSIF  D(U,.K. 

GO  TO  190  BeStA 

ISO  DO  160  1=1, hL 

READ  (5,7100)  Hb ( I ) , RHOM ( 1) 

REAO  (5,7100)  EX(I),ET(I)»XXNU(I)» THNU (I ) » GX T £  I) 
160  REAO  (5,7100)  SaTCI).SAC(I) 

IF  (KTYPF.NE. 3. AND. KTYPE.NE.fi)  GO  TO  190 
IF  (KCAP.Nfc .0)  GO  TO  190 
REAO  (5,7100)  EC , GC , OC 
190  REAO  (5,7000)  fJSPR.  NTSC.  (*F  I  (I ) ,  1  =  1 ,  a) 


Reproduced  From 
Best  Available  Copy 


IF  (NSPR.EG.O)  GU  TO  205 
DO  200  1=1 , NSPR 

200  READ  (5.7110)  10 IR ( I ) . NSPG (I) . NSPB ( I ) , BIGK  (  I ) 
205  IF  (NTSC. EG. 0)  GU  TO  211 
00  210  1=1, NTSC 

210  READ  (5,7120)  iSLUC ( I ) . K TSGB ( I) , B I GC ( I ) 

211  NIPF  =  0 

00  213  1=1,4 
NFIX  =  NF 1(1) 

IF  (NFIX. EG. 0)  GO  TO  213 
NIPF  =  NIPF  +  1 
NFIXX  =  IaBS(NFIX) 

REAO  (5,7100)  (TIP(J,I).J=1,NFIXX) 


REAO  (5,7100)  (FIP(J. I) , J  =  1 FIXX) 
NSF  =  N«AR 

IF  (I.GT.2)  NSF  =  MBA  W 
NSFX  =  NSF 

IF  (MHx.LT.C)  NSFX  =  1 
READ  (5,7100)  (SF1P(J,I),J=1 .NSFX) 

IF  (NSFX.GT.l)  GO  TO  213 
UO  212  J  =  2 , N S F 

212  SFIP(J.I)  =  SFiP(l.I) 

213  CONTINUE 

KSTIF  =  IABS(NSG)  ♦  I AftS (VSB) 

IF  (KSTIF. EG. 0)  NCOUP  =  0 
LS6AK  z  1 

IF  (NSHEAR.FG.OJ  GU  TO  216 
REAP  (5,7100)  C.GG.&R 
REAO  (5,7000)  lkBak.lSBaP.LPIJT 
REAO  (5.7000)  (LRG(I } , 1=1 .LRfcAR) 
READ  (5. 7000  (LSBCI)  ,Isl.LS6AP) 

IF  (LCDT.EG.O)  GU  TU  216 
00  21  a  Tzl.LOUr 

210  READ  ( 5 »  7  0  0  0  )  L Ml/U T  ( I  )  ,  L SOP T  f  T  ) 

216  bUCKS=0.0 

IF  (NBl'CK  .GT.  0)  KEAO(9,7!00)  «L'C*S 
IF  (KSTIF. EG. 0)  GU  TO  000 


IF  (MSG. EG. 0)  GU  lu  260 
\SGY  =  IAHS(nSG) 


READ  (5,7000)  («SG( I ) , T=1 , mSGY) 

NSGX=NSG 

IF  (wSG.LT.O)  uSGx  =  1 
OP  220  1=1, NSG* 

READ (5, 7  1  no )  SI  DEG ( I ) , ESTPG ( I ) . GR ARG ( I ) , RH0S TG  ( 1 )  , 
1  SIGubT ( I) .STGOGC (  I) 
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ETSTRG(I)  =  0. 

EPSG ( I )  =0. 

IF  (K0ERV.EU.2)  READ  (*5*7100)  ETSTRG ( I) , EPSG {  I) 

REAO(5*  7000)  MSEGGd  )  .KSUPGtl) 

IF  (SIOEGfl)  .Ed. 2.0)  KSUPGCI)  =  0 
NS  s  MSEGG(l) 

READ  (5,7100)  SIGJG(I) ,HOG(I) , AwEeG(I) 

DO  220  1=1, NS 

RE AO  (5,7100)  HSTGCL. I).BSTGfL.I) 

220  CONTINUE 

IF  (NSG.GT.-2)  GO  TO  260 
00  250  I=?,NSGY 
SIOEG(I)=Sl OEGll) 

ESTRGCI)=6STRGin 

GBARG(I)=G8ARG(1) 

RH0STG(n=RH0SlG(1  ) 

SIGOGT  ( 1 1  =  SlGuGfd) 

SIGOGC(I)  =  SIGOGC(l) 

ETSTRG(I)=£rSTWG(l ) 

EPSG  ( I )  =  EPSG ( 1 ) 

NSEGG ( I ) snSEGG ( l ) 

KSUPG(I)  =  KS'jPGll) 

NS=NSEGG(  1  ) 

BIGJGt  T)=8IGJG(1) 
hOG(I)=HOG(l) 

00  ?a  0  L  =  1 ,NS 
HSTG ( L , I)s“STG(L, 1) 

200  8STG(L, T1=! ;STG(L,  1  ) 

250.  CONTINUE 

260  IF  (NS8.FQ.0)  GO  TC  ann 
NSrtr  =  IARS(N$3) 

KfcAO  (5,7000)  (ft Sd ( I ) , Ts 1 , NS8  Y ) 
nSh  X=NS* 

IF  (nSR.LT.  0)  "•SOX  =  1 
00  300  I  =  t,r*S3X 

RE  A  0(5, 7  1  OO)  SI0fce(n,FSTR3m.GRAh8(I),RM0STB(I), 

1  SIGU6 T ( I ) . SIG06C (I) 

ETSTHB(I)  =  0. 

EPSB(I)  =  0. 

IF  (NOE  W.E0.2)  REfvo  (5,7100)  ETS  TPH  f  I )  ,  EPS6  ( I ) 

KE A 0(5, 7000)  fJSEGn  (1  )  ,KSUPB  (  T)  ,  NSTH  ( I) 

IF  (S1DERC I ) .EN.2.0)  KSUPfi(I)  =  0 
NS  =  NSEGB(I) 
nSTBXsnSTR  f  I ) 

00  30  P  K  =  1,r;STyx 

READ  (5,7100)  oIGJUR,  I  )  ,HC>R  (K  .  I)  ,t«tTC  (K  ,  I)  ,  A/^ERBCK,  !) 

00  30<>  l.=  l,NS 

RE  AO (5 , 7100)  H5T6(L,K.I),*jSTR(L,K.n 
300  CONTINUE 

IF  (fSH.GT.-2)  GU  TO  ann 
00  350  I=2,';5ht 
S I  OF  B (  1)=S10Eh(1  j 
ESTRR ( r ) =E ST«^(1 } 

6  a  A  ft  F>  (  I  )  =  P  h  A  ^  H  (  1  ) 
ft^OSTB  {  n  =  ft  H  0  S  T “ ( 1  ) 


From 

Best  Available  Copy 
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3IG09T  (I)  =  SIGOGT(l) 
SIGUeCCI)  *  SIGOBC(l) 
ETSTRStnseTSTHttU) 
EPS8(I)  s  EPSB(l) 

NSEGB ( I ) SNSEG8 ( 1 ) 
KSUPB(I)  =  KSUHB(l) 

N3T0  Cl ) sNS  T0  ( 1 ) 
NS=NSEG8M  ) 

NST8X  =  NST8(t) 

DO  330  K=l,N3TbX 
0IGJ6(K,n=8lGJa(K,l) 
H08(K,I)=H08(K,  1) 
«ETC(K, I)=9ETCl«,  1) 

00  330  L=1»«S 
H3TB(L,K, 1)SH3T6(L,K.  1) 
330  RST0(L.K.I)=8ST6(L,<.1) 
350  CONTINUE 


Reproduced  From 
Best  Available  Copy 


400  READ  (5.7100)  ( t FG C I . J ) . 1  =  1 . MB ) . J  =  1 , MG ) 

READ (5#  7 100)  0ELT1*,  TSTOP.  PRINT 
IF(IN'OUT.EG.O)  GO  TO  2100 
PRINT  OUT  THE  INPUT 
«RITE(b»7l70) 

WRITE  (6,7200)  ^G,M8.MRAR.NBAR,L8Aft 
*RITE  16.7210)  tMGM(I).I=l,PG) 

WRITE  (6.7220)  (*&*( l) • 1st ,M8) 

WHITE  (6. 7  225)  i'«3YvG»NSYNB»NPLTeNBUV,!Mt  '<D»NDERV 
WRITE  (6.7600)  NSG.NS6.KCOUP 
WRITE  (6.7150)  «uOUT 

IF  (NNOUT.GT.O)  wRI  TE(6.7180)  (POUT  (I )  ,  NQUT  ( I)  .  1  =  1  ,  NiMQUT) 
WRITE  (6.7185)  NKp 

IF  (NKP.GT.O)  -RUE  (6.7180)  (  KPG  ( I )  .  KP8  (  I )  .  I  =  1 ,  NK  P  ) 

WRITE (6. 1 1600)  NL, NSHEAR.KBUCK . XIP 
IFCNPLI.EO.O)  v RI fE(6.7?3C)  THE! AO 
IF(NPLT.ER.l)  .RITE(6.7260)  THE! AO, A 
IF  (N0ERV.EQ.2J  GO  TO  1180 
00  1160  1=1, NL 

1160  WRITE  (6.  1  1700)  I . HM ( I ), RHQM f T ). E X ( I ). ET ( I ). X XNU C I ), T HMJ  ( I )  , 
1  GXT(I),SAT(I),SaC(I) 

IF  (KTrPE.NE.3.AND.<TYPE.6E.«)  GO  TO  1190 
IF  (KCAw.nE.OI  GU  TO  1190 
wRITE  16,11900)  EC.5C.0C 
GO  TO  11 Q0 

1180  wRITE  (6.72«0)  (h«- (I  )  ,  rhUP  (  I)  .  E«  ( I )  ,  1  =  1  #  NL) 

WRITE  (6.7300)  Tno.SIGO.EP.EPSTF 
1190  WRITE  16.  12000)  "-SPR 

IF  (NSPR.EU.O)  GO  TO  1210 
WRITE  (6.12200) 

00  1200  I  =  l,.'iSPR 

1200  WRITE  (6.12100)  IDIR(I)  , NSPG(I)  , NSP0 (  I)  .8IGK  (I) 

1210  WRITE  (*,12300)  uTSC 

-RITE  (6.12600)  (f FI ( I) . 1=1 ,4) 

IF  (NTSC. EG. 0)  GU  TO  1225 

wRITE  (6,12400) 

I'O  122  0  1  =  1  , F SC 

122"  ^RITE  (6.12500)  I SLUC (I) . K TSGP ( I ) . 8  I GC  ( I ) 


115 


1225  IF  (nipf.EQ.O)  6u  TO  1240  Reproduced  From 

DO  1230  1  =  1,4  Best  Available  Copy 

NFIX  =  I4RSCMFHI)) 

WRITE  (6,12700)  It  (TIPCJ. I). Jsl.NFIX) 

WRITE  (6,12800)  (FIP ( J , I) , J=1 , WF IX ) 

NSFX  =  N9AN 

IF  (I.GT.P)  NSFX  =  MBAR 

IF  (NFI (I) .LT.O)  NSFX  =  1 

WRITE  (6,12900)  (SFIP ( J, I ) , J=1 , NSFX  ) 

NFI(I)  =  NFIX 
IF  (MFIX.GT.l)  GU  TO  1230 
TIP(2, I)  =  1.E6 
FIP (2, I )  =  FIP(1,I) 

1230  CONTINUE 

1240  IF  (MSHFAR.tM.O)  GO  TO  1245 
wRITE  (6,13000)  C,GG,G0 
WRITE  (6,13100)  LHSAR,LS0AR,LOUT 
WRITE  (6,13200)  (LRG ( I) , 1=1 » LR8 AR ) 

WRITE  (6,13300)  (LSfl(I),I=l.LSRAW) 

IF  (LOUT.GT.n)  WRITE  (6,13400)  (L ROU T ( I ) , LSOUT ( I ) . I = 1 , LC U T ) 
1245  IF(NBUCK  .GT.  0)  ARITEC6, 13500)  BUCKS 
IF (NSG  .EG.  0)  GU  TO  1300 
WRITE  (6,7700)  (KSG(I) , 1=1 ,NSGY) 

DO  1250  Is l f NSG X 

WRITE (6,7800)  SlOEGd ) , ESTRG  f I ) ,GBARG( I) , RHOSTG ( I) , S I GOG T  ( 1)  , 
1  SIGOGC(I) 

IF  (NPERV.EC.2J  WRITE  (6,7830)  ETSTRG ( I ) . EPSG (I) 

WRITE  (6.7  850)  f-SEGG(I)  ,KS<JPG(  I)  .BIGJGCI )  ,H0G(1 )  ,  Av  EwG  (I) 

MS  =  NSEGG(I) 

DO  1250  1=1, NS 

WRITE  (6.7900)  MSTG(L, I),eSTG(L,I) 

1250  CONTINUE 

1300  IF  (MSB.EQ.O)  GO  TU  1400 

WRITE  (b.  80(10)  (KSb(I)  ,1  =  1  ,i*88Y) 

Dn  1350  isl.NSBX 

WRITE (6, 8050)  SlUtB(I) .ESTHBCI) ,GBARB(I) ,RHOSTB(I) , SJGOPT(I) , 
1  SIGOBC(T) 

IF  (MDERV.EW.2J  WRITE  (6,8070)  FTSTRU (I) , EPSB (I ) 

WRITE  (6.8090)  NSEGB( I ) .KSUPB(I) ,NSTB ( 1) 

NS  =  MSFGB(l) 

NSTBX=NSTB  (I) 

DO  1350  K s  1 ,  f-. S  T b x 

WRITE  (6,8300)  K,bIG J8 (K , I ) .HOB (K, T) ,HETC (K, I ) , A WE BE (K , T ) 

DO  1350  L=l, MS 

WRITE (6, 7900)  HSTBCL.K.T),9ST9(L,K,n 
1350  CONTINUE 

1400  WRITE  (6,7400)  (  ( F  G  (  I  .  J )  .  1=  1  .  <»8  )  ,  J=  1  .  MG) 

WRITE (6. 8200)  UtL T 1* , TSTOP , PW  PiT 


2100  1=0 

MGNRsO 

DO  2150  M=l,i’G 
R-.V  s 

DU  2150  N=t,"R 
NM  =  MSN  CO 
MUSE  (N,  0  =  1 
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Reproduced  From 

IF(I.EQ.nnout)  go  to  2130  Best  Available  Copy 

00  5110  J  =  1  # NNOUT 

IF (MW. EG .MGUT ( J ) • AND . NN.EG . NOUT  CJ)  )  GO  TO  2150 
2110  CONTINUE 

GO  TO  2130 
2150  MUSE (N » M) sO 
1  =  1  +  1 

GO  TO  5150 
2130  MGMB=MGMR+1 
2150  CONTINUE 

IF  (NSHEAR.EO.O)  GO  TO  2500 
I  =  0 

DO  2400  Msl.LRBAR 
MM  =  LHG(M) 

00  2400  Msl,LS8AR 
NN  =  LSB(N) 

LUSt (N#  M)  =  1 

IF  U. EG. LOUT)  GO  TO  2400 

DO  2200  J= 1 »  LOUT 

IF  (MM.6Q.LR0UTCJ) .ANO.NN.EQ.LSOUT(J) )  GO  TU  2300 
2200  CONTINUE 

GO  TO  2400 
2300  L'JSE(N.M)  =  0 
1=1  +  1 
2400  CONTINUE 
2500  CONTTNIJF 

mgmB2=2*mghb 

XJS180.O/THETAO 

IF(NPLT.EO.O)  XJ=PI/THETAO 

NBMDt  s  MBND/10 

N0NO2  5  NRND  -  1O*N0NO1 

NBUVl  =  N8U0/10 

N3UV2  =  NBUV  -  1 0*NBUV 1 

RETURN 

7000  FURMA  T (6112) 

7100  FORMAT (*F 1 2. 1 ) 

7110  FORMAT  (3U2.F12.1) 

7120  FORMAT  12112»F12.1) 

7150  FORMAT  (BHONNOUT  =  133 

7170  FORMAT (lObMl INPUT  DATA  FOR  OEPRCP  (MODIFIED  TO  INCLUDE  EXTRA  8.C., 

i  el  a  st  r  c  a  no  tursiunal  springs  and  in-plane  forces)) 

7100  FORMAT  (514) 

7105  FORMAT  (7H0-VKP  =  1 3 » 

7200  FORMAT (  10H0MG  =  I5/IDH  MR  =  I 

12/10H  mHAR  =  I2/10H  NBA R  =  I2/10H  L8AR  =  15) 

7510  FORMAT  (10HOMGN  =  (1015)) 

7520  FORMAT  (1OHON0W  =  (1015)) 

7555  FURmaT  (IOhonSVmG  =  I5/10H  N  S  Y  M  8  =  I2/10HONPLT  =  15/ 

5  10m  NBUV  =  12/ 

3  10H  NHN[)  =  15/  10H  NDERV  =  15) 

7530  FORMAT (1 7H  THEIao,  [m  =  Elb.fl) 

7260  FORMAT  ( 1  7  m  THE  T  AC  r  OEG  =  E16.B/17R  A,  IN  s  E16.B) 

7240  FORMAT  M2HP  hm.  [m,  4*  ,  51  hR'HQM,  LP-SEC**2  / 1  N**4 . 4  x  , 

1  7hEm.  PS?/(3E17.3) ) 

7300  FORMAT (17H0T  Hi  =  E16.B/17W  SlGOf  PSI  =  E16.P/17H  £P, 
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Reproduced  From 
Best  Available  Copy 


1  PSI  =  E16.8/17H  EPSIF.  IN/IN  =  Elfa.8) 

7400  FORMAT (feHOFG  =  /(5E14.6)) 

7600  FORMAT  (10HCNSG  s  I2/10H  NSR  =  I2/10H  KCOUP  =  T2) 

7700  FORMAT  (tOHOKSG  s  (1015)) 

7800  FORMAT  (23H05 IDES  =  FU. 1/ 

2  23H  ESTRG,  L0/IN**2  =  F15.6/23N  GBARG.  LB/IN**2 

3  23H  RHIJSTG#  LB/SEC2/IN4  =  E15.6/23H  SIGOGT,  LB/lN**2 

4  23H  SIGOGC,  Lb/lN**2  =  E15.6) 

7830  FORMAT  (23H0ET8TRG,  L8/IN**2  =  E15.6/23H  EPSG.  In/IM 

ieis.6) 

7850  FORMAT  (ION  NSEGG  =  I3/JOH  KSUPG  =  13/ 

1  23H  BIGJG,  IN**4,  s  E15.6/2 ?H  HOG,  IN 

2  22HA>vFHG,  In**2  =  E16.6/ 

2  15X,9H  HSTG,  In,21X,8HBSTG.  IN) 

7900  format  (2  (tax  ,ElS.e>n 

8000  FORMAT  (10H0KSB  =  (1015)) 

8050  FORMAT  (?3H0SlL)Ed  =  F4.1/ 

2  23H  ESTRB,  L'3/1N**2  =  F15.6/23H  G8ARB,  LB/IN**? 

3  23H  RHOSTR,  LB/SEC2/IN4  =  E15.6/23H  SIGOBT,  LB/IN**2 

4  23R  SIGOBC,  LB/lN**2  =  E15.6) 

8070  FORMA)  (23H0ETSTRH,  (_8/IN**2  =  E15.6/23H  EPSR,  IN/IN 

) E 15.6) 

8090  FORMAT  (10H  NSEGB  =  I3/10H  KSUPH  =  I3/10H  NS T B  =  13) 

8200  FORMAT (15H00ELTIM,  SEC  =  E16.8/15H  TSTOP,  SEC  =  E1&.8/15H 


1  =  E16.8) 

8300  FORMAT  (9H0ST  A  T  i(JN  1 2 /5 X , 1 5HB I G JR ,  If  **4  =  £15.6/ 

1  5X,14HM0B,  In  s  E15.6/5X, 14RBETC,  IN  s  E15.6/ 

2  5X  ,  1  un  Ak-EBB,  1 1  j  *  *  2  E 1 5 . 6/ 

2  15X,9h  wSTR,  IN,21X.8HBST8,  IN) 


1600  FORMAT  MOHO-'iL 


=  12/1 OHOWSHE  Aft  =  12/ 


1 


10HON8UCK  =  I2/17H0XLP,  I M 


=  E  1 6 . 8 ) 


1700 

FORMAT 

(6H0LA  rtRI3/27H 

Hfl  , 

IN 

1 

27H 

RHGM,  L6-SEC**2/If **a 

=  £16.8/ 

2 

27H 

EX,  PSI 

=  £16.8/ 

3 

27H 

FT,  PSI 

=  Fib. 8/ 

4 

2  7  H 

XXNlj 

=  £16.8/ 

5 

27h 

T  H  U 

=  £16.8/ 

6 

27H 

G  X  T ,  Pol 

=  £16.8/ 

7 

27H 

SAT,  PSI 

=  £16.8/ 

8 

27H 

SAC,  PSI 

=  £16.8) 

=  E 1 6 . 8/ 


1900  FORMAT  MIHOEC,  RSI  =  E16.S/11H  GC,  PSI 
1  11H  DC,  IN  =  E  1*?.9) 


=  E16.8/ 


2000  FORMAT 
2100  FORMAT 
2200  FORMAT 
2300  FORMAT 
2400  FORMAT 
2500  FORMAT 
2600  FORMAT 
2700  FOR* AT 
2800  FORMAT 
2900  FORmaT 
3000  FORMAT 
3100  FORMAT 
3200  format 


fBHOrj.SRR  s  131 
<5X,3Io,E15.6) 

( 3  7  h  I  r>  T  R  NSPG  NSPB  K  (LB/IN)) 

(8 HO NTSC  =  (3) 

( 1 1  X,/jhCU(.'E,2X,8HL0CA  TI0f:,  3>  ,  1  iHC  C  IN-LB  /RAO)  ) 
(5X,2I6,E15.6) 

(HHOMFl  =  414) 

(^nOrifi-jKOaur  I3/12H  TJMFS  =  5E1  3.4/H  «  *,5F  1  3.4)  ) 

(lcti  FORCES  =  5E  1  3.4/ (  l  1  X  ,5F  l  3.U)  ) 

MQM  SCALE  FACTORS  =  5E13.4/(18X,5E13.4) ) 
f  8H0C  =  E  15. 6/6 H  GG  =  E15.6/BN  GB  =  FIS.r) 

(  1  G  R  0  L  *  t>  A  K  =  T2/JOH0LSHAR  =  I2/10H0LUUT  =  T2) 
rwHOLKb  =  (1015)) 


3800  FfiRM^T  (*h  LSh  =  (  )  0  T  5 )  ) 


E 1 5. 6/ 
E 1 5. 6/ 


E  1 6 . 6/ 


E 1 5. 6/ 

E15.6/ 


PRINT 
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4 


4 


3400  FORMAT  ( ifcHOLKUUr,  LSOUT  =  /(5X.2I5)) 
3500  FORMAT (10H08UCKS  =  615. fei 
END  ' 

H  s  HM  (Ml.) 


F3sRH0/CJH3 

F1sOM22*A**3/F3 

F2=F/F3 

F  4  =  F3/fA*CMU) 

F5  =  OMi i*A**3/F3 
F3=F3/CA*CM22) 

DT  =  OELTIM 

CALL  OTSTEP  (FI , F2 » F 3 » F 4 » F 5) 
IF  COT. GT. 0.0)  OELTIM  s  DT 
NELP  =  1 


Reproduced  From 
Best  Available  Copy 


NZP  s  2 

IF  (KTYPt.E0.5J  GO  TO  2730 

NLZ(l)  =  1 

ZC(1)  =  -HBAR 

MLZ(2)  =  t 

ZC  (2)  =  zcm  +  H 

IF  (KTYPE.LT. 3)  GO  TO  2745 

NLZ (2)  =  3 

ZC(1)  =  ZC(1)  ♦  .5*hm( 1 ) 

ZC  (2 )  =  ZC(2)  -  ,5*(HM3)  -  HM(2)) 
GO  TO  274S 
2730  HT  =  -HHA« 

HS  s  HM(l) 

00  2740  1=1, ML 
ML Z ( 2* I  -  1)  =  1 
NLZ(2*IJ  =  I 
ZC(P*I  -  1)  =  HT 

IF  (I.GT.l)  hS  =  H*(I)  -  MMCI-l) 

HT  —  HT  +  HS 


2740  ZC l 2*1 )  =  HT 
NZP  =  2*N(, 

2745  IF  (KDAM.E0.2)  GO  TO  2800 
00  2747  T=1,NL 
TCRIT(I)  =  SAT  1 1 ) 

2747  CCRIT(T)  =  SAC(I) 

2750  CONTINUE 
2«o0  HETURn 
3000  K£RR  s  i 
RETURN 
END 

* 0 E C K  OSET3 

SUBROUTINE  DS£  f  3 


*C ALL  C8LK1 
♦CALL  CBLK2 
♦CALL  CHLK3 
♦CALL  CBLK4 
♦CALL  C8LK5 
♦CALL  CHL«7 

♦  CALL  CBt.KP 
♦CALL  CHLK10 

♦  CALL  C  8  L  K  1  i 
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Reproduced  From 
Best  Available  Copy 


♦CALL  C8LK13 
•CALL  CBLX15 
•CALL  CNOVA 
•CALL  C8LANK 
•CALL  C8LKJ8 

PRINTOUT  DESCRIPTION  OF  OEPRL'P  DATA 

2800  fcRITE (6,9300) 

IF(NPLT.EG.O)  WK>IT£{6, 9400) 

IF(NPLT.EQ.l)  wk1TE(6,950(.) 

GO  TO  (282 0,284 0,2860, 2080,2900) ,  KTYPE 
2820  WRITE  (6.9600) 

GO  TO  2950 
2840  WRITE  (6,9700) 

GO  TO  2950 
2860  WRITE  (6,9800) 

GO  TO  2950 
2680  WRITE  (6,9820) 

GO  TO  2950 
2900  WRITE  (6,9840) 

2950  4RITF  (6,13500) 

IF  (NBUVl.EO.l)  wRHE  (6,  13600) 

IF  (NBUV1.EP.3)  WRITE  (6,13700) 

IF  (NRlJVl  .FG.5)  WRITE  (6,1  3800) 

WRITF  (  6 » 13900J 

IF  (NBUV2.6U. 1 )  WRITE  (6,13600) 

IF  (MHOV2.EfJ.3J  .RITE  (6,13700) 

IF  (MBUV2.EG.8 j  WRITE  (6,13800) 

-RITE  (6,14000) 

IF  (GBND! .EG. 1 )  WRITE  (6, 9900) 

IF  (NHNPt.EQ.2)  WRITE  (6,9920) 

IF  (bBNI)l  .E(i. 3)  WRITE  (6, 993") 

IF  (NHMP1.EG.4J  WRITE  (6,9940) 

IF  (NBND1.EG.5)  .'RITE  (6,9950) 

IF  (NPMC1  .GT.5)  wRUE  (6,9955) 

IF  (M0MO2.EG.1)  WRITE  (6,9960) 

IF  (NBND2.EQ.2)  write  (6,998") 

IF  (NPND2.EQ. 3 )  WRITE  (6,9990) 

IF  (MHN02.EG.4J  « R I T  E  (6,1  0000) 

IF  ( i  J  B  f  i  0  ?  .  E  G  .  5  )  wrUE  (6,  10050) 

IFCMBM02.GT.5)  WRITE (6,10055) 

I F  ( NDER  V  .  E  (•’  .  1  )  ww  l  T£  (  6 , 1  0 1  0  0  ) 

IF (MOFRv .EG. 2)  wk I Tt (6, 10200) 

IF  (KCOUP.EG.l)  write  (6,102501 
WRITE  (6, 10800)  f-’G, -*»»), mbaR.MHaK.LBAR 
WRITE  (6,1 0820 J 
DO  2970  *=1,VG 
MH  s  MfJM(H) 

DO  297  0  M  =  1  ,  mh 
MM  =  fJpMfiv') 

IF  («USF (M,M) ,b9.0 )  GO  TO  2° 7  0 
-RITE  (6,  IGB^OJ 
2970  CONTINUE 

IF  ( MS HF  A R  , F .10 .  u  .  n i v  f).  (KTYPF+1  ) /2.EK.2)  WRITE  (6,14100) 
IF  (("Shear. nf  . oj  -  riie  (6, 14200) 
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WRITE  (6,1 0850)  XLP 
TF(NPLT.EG.O)  wRI IE(6» 10900)  THETAO 
IF  (NPLT.EO*,  1 )  w«ITE(b,l!OrO)  THETAO,  A 
IF  (N0ERV.EO.2)  60  TO  3500 
WRITE  (6,12050)  H8AR. (I, Isl,NL) 

WRITE  (6,12100)  (HM(i),isl.hL) 

WRITE  (6.12200)  (KH0M(I).1=1.NU  RCPT°  -Vab\C  COpV 

WRITE  (6,  12300)  (EX  ( I ) ,  1  =  1 » ML)  Best*va 

WRITE  (6.12000)  (F.T(I)  ,1  =  1  ,fJL) 

WRITE  (6.  12500)  (XXN(j(I),Isl.NL) 

WRITE  (6.12600)  ( THNU ( I) . 1* 1 . ML ) 

WRITE  (6,  12650)  (6X T  (  I)  .  T«1 . «U 
IF  (KTYPF.ME.t .ANU.KTYPF.MF.3)  60  Tu  3300 
WRITE  (6,12000)  (3AT( I) . 1=1 »NL) 
wHITF.  (6.  13000)  (SACm.Ul.MU 
GO  TO  3«0O 

3300  WRITE  (6.12700)  (3A T ( T) . T=1 . NL) 

WRITE  (6.12800)  ISAC(I).I=1.NL) 

3a00  IF  (K1  YPE.ME.3.AN0.KTYPF_.NE.<O  GO  TO  3600 
IF  (KDAM.NE.O)  GO  TO  3600 
WRITE  (6.13100)  EC.GC.OC 
GO  TO  3600 

3500  IF  (KTYPE.E0.3)  wH I TE  (6.13300)  HPAP 

WRITE  (6.11100)  h.RHO.EL.TMU.SIGO.EP.EPSIF 
WRITE  (fc. 13400)  NSPR , NT  SC 
3600  WRITE  (6,1  1200)  UF5(a..M),N=1,HR),M=1,MG) 

WRITE (6, 1 1300)  0ELT1M.TST0P. PRINT 
IF  (MPERV.EO.l)  60  TO  4020 
00  4000  Ksl.LHAR 
ZH(K)sGX(K ) 

ZF (K)sZH( K) /A 
ZG(K)=GX(>0**2 
a000  CONTINUE 

Z  A  {  1  )  =  Z  H  (  1  )  /  A 
ZA(2)  =  ZRC2)/A 
INZ(I)  =  l 
INZ(2)  =  I  BAR 
4020  NMSYMG  =  0 
MNSYNH  5  o 

IF  (NPN01 .GT .3)  wMSYMG  =  1 

IF  (NHM02.GT.3J  KNSr*B  =  ! 

IF(NHN01 .EO.b)  WnSYmGsO 
IF (MBNC2.E0 .6)  wwSYMBsO 
IF  (MHUV1.GT.3)  MnSYmG  =  1 

IF  (NHl.iV2.GT. 3)  nnSTT^B  s  1 

IF  (NNSYMG.EM.  1  .AWU.t'JSYMG.EG.O)  *PITE  (6,13200) 

IF  (MNSYWB.EN.l.A-'.-O.NSYMB.fcu.Q)  WRITE  (6.  1320C) 

NGT  s  mb ar 
NBT  =  NEAR 
NG  =  fGT 
MB  =  MB  T 

IF  (WSYMG.EU.l)  MG  =  (NGT+1 ) /? 

IF  (NSYmh.Eu.  1  )  NH  =  (NBT+U/2 
M  Y  j?  =  3*mgmb 

HIM  =  P  I  /F  LG  A  T  ( e!*  LMB  AR®  1  )  ) 

PIN  =  PI/FLOA T (2*  (MBAR-1 ) ) 
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IF  (NSYNG.EQ.l)  PI'*  =  2.*PJM 
IF  (NSYM8.eQ.lt)  PIN  s  ?.*PIN 
RsA/H 

BUCK  =  «ABS (BUCKS! 

BUCK9s-ABS(8UCKS)/R 
XL-XLP / ( P I  *  A ) 

XL  1*1 • 0/XL 
XL2*XL**2 
XL3*2. 0*XL  1 
XL4*2.0*XL2 
XL5=XL«*R 
XL7S1.0/XL2 
CN8  =  XL  **2 

IF  (KTYPE.tU.3.  Ai40.N0ERV.EQ. 2)  CN8  =  2e0*CN8 
CN9  =  CN8/(2.0*K) 

IF  (NDEBV.EQ.l)  GO  TO  4040 
CN10  =  2.0*CMft 
CN11  =  CM9/(3.0*K) 

GO  TO  4050 

404C  CM 1 0  =  2.*CM»*« 

CM  11  =  CM  10 
4050  XJ2  =  X J**2 
XJ3sXJ*XLl 
XJ4=2.0*XJ 
xJ5=2.o*x J3 
OPRT  1=PKIiiT*DElT1'-' 

XLP1=XL5 
XLP2=2.0*XL*R 
IF  (M&PW.fcU.O)  GU  TO  4075 
XLP  3  =  2.0*XJ/(A*t-’I*H*XLP) 

IF  U'SYKT,*f:SYfb.tU.n  GO  1C  4075 
ACCOUMT  FUR  SYMMETRY. 

XLP 3  =  4.0*XLP3/  FLOAT  (  (N?ynr,+  1  )*  (MSYMB+1  )  } 

00  4060  L=lfMSPK 
I  s  MSPGfL) 

J  =  MSPPfL) 

IF ( I.EO.mBAK)  nlGK(L)  =  0 . 5*R TOK (L ) *FLO A T (NS YMG+ 1) 

IF  (J.EO.rHAk)  nlGK(L)  =  <<  .  5*H  T  GM  L  )  *F  LO  A  T  (  NS  Y  Mb ♦  1 ) 

40b0  CONTINUE 

4075  IF  (NIPF.EO.O)  GO  TO  4155 

XLPfe  =  2.P*(PI*xL)**2*B/(A*THtTA0*XLP**3) 

IF  (UPLT  .t'l.  1)  XLP  6  =  XLPb*140.0/PT 
00  41S?  1*1,4 

IF  (OF  I  1 1)  .fc'4.0)  GO  TU  4152 

IF  (  I.LT,3.AiviO.NSYMB.EQ.O)  SFIP(NBAR,TJ  =  0.5*SF  TP  (MfjAR  ,  I  ) 
IF  U  .GT.2.  AwD.MSY^G.EO.O)  SF  T  P  (MBAP  ,  I  )  =  0 . 5*  SF  J  P  (  MB  A  R  ,  I  1 
JLIPU)  =  2 
4  152  CUM  If  HE 
4155  CONTimiF 

SIMPSf'M  s  KHLF. 

MQ A H  ai  d  f.HAb  ■«US  T  BE  ODD  NUMBERS  FUR  FULL  PAuEL. 
on  4  1  EC  T  =  l,klBAR 
F  -  [  *■  1 

■  GAM(I)  =  F  *  P  I 

kGU)  =  GAUUJ*XLP/PI 


iced  From 

Reproducer 

BestAvailaWe 


* 


."f 
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H J ( I )  =  2.0 

pima (i)  =  p„0*pim/ (3.o*xlp**2)  Reproduced  From 

IF  CNSY«G.EG  1)  PIMACI)  =  .5*PIMA(I)  Best  Available  Copy 

IF  (CI+1J/2.EQ.1/2)  HJ(I)  =  0.0 
4160  CONTINUE 

HJ(1)  =  t.O 

IF  (NSYMG.EQ.l)  HJ(mqar)  =  1.0 
IF  (NSYMG.EG.OJ  PlMA(MBAR)  =  PlMA(MBAR)/2.0 
00  4200  Isl.KiBAP 
F  =  I-S 

BETRC1)  =  F*PIN 

XB  (  I  )  =  RETR(n*TH£TA0/Pl 

HK II)  =  2.0 

PINA(l)  =  2.0*PIN/3.0 

IF  (NSYMB.E0.1  )  PltMA(I)  =  .5+PINA(I) 

IF  (  n  +  1) /2.fcfJ.l/2J  HK  Cl)  =  4.0 
GAMHA(I)  =  XB ( 1 )  -  .5*THETA0 

IF  (NPLT.E0.1J  GAMMA (I)  =  GAMMA(I)*KFR*PI/180. 

4200  CONTINUE 

IF  (MPLr.FU.l)  1HETA0  =  THE T A  0 * P I / 1 H 0 . 

NMASS  =  NBAft 

mm ass  =  r bar 

HK (1)  r  1.0 

IF  (NSYmh.EG.1)  hk(nBAR)  =  1.0 
IF  (fJSYMM.EG.OJ  F*  1 1>«  A  ( NH A  ft )  =  PIFA  f  NH  ARJ /2. 0 
DO  4 2 SO  1=1, M RAP 
4250  PIMA(I)  =  PIHA(1)*MJ  (I) 

00  0300  1  =  1  ,  PR  A  k 
4300  PIMACI)  =  P INA ( 1 ) *HK ( I ) 

DO  4400  I=l,M0Ak 
00  4400  J  = 1 ,  N  R  A  ft 
44  00  NUSE  (J,  I)  =  2 
II  =  O 

00  4430  r  =  l,'‘GT 
DO  4470  J=l,f.r>l 
IF  (11 .FQ.NKP)  GO  TO  4430 
00  441  0  K  =  l,r»KP 

IF  (  I.EO.KPGIK  J  .  At,0.  J.ECJ.f  PB(K)  1  GO  To  4«2(t 
4410  CONTINUE 

GO  TO  4430 
4420  NUSt (J,  I)  =  3 
11=11+1 
4430  CONTINUE 

CALL  ST  IFF  (  1  ) 


IF  (ivTSC.EU.O)  bU  TO  4700 
TORSIONAL  SPRINGS. 

XLP4=FI**?/(3.*n*FL0ATfN8AP-t}*xLP**3) 

IF (NSVMH.Fn. n  XLP4=2 . 0*XLP4 

XLP5=f  P  I  *X  J*XLJ  *-*2  /  (3.»H*FL0ATf  WRAH-I)  *A*THETAO*XLP  +  *?J 
IF (NSymg.EG.1 )  XLP5=2.0**LP5 
TGI  =  0. 

TG?  = 
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4060 

4063 

4065 


u  0  88 

<,090 

4  0  42 
4094 
4096 

41  CO 

•4110 

4115 

412  0 

4125 

4  1  X  u 
4  1  <5 


TB1  =  0. 

T82  =  0.  Reproduced  From 

n  T  sc  x  =  o  Best  Available  Copy 

DO  4065  T  =  1 ,NTSC 
J  s  1 

K  a  ISLOCCl) 

IF  (CK  +  D/2)  4060,4083.4065 
J  =  MHAR 
GO  TO  4085 
J  =  NR AR 

NTSCX  =  NTSCX  +  J 
J  =  NTSCX  ♦  1 
DO  4096  IT  =  1 , NT  SC 
I  =  NTSC  -  II  ♦  1 
K  a  ISLOC(I) 

KK  a  1 

IF  (K.LT.O)  KK  =  ‘*8  AR 
IF  (K.LT.-2)  ^  =  MRAR 
IF  CK.GT.O)  KTSGBCJ-1)  =  KTSGR(I) 

I AK  =  IAHS(K) 

DU  0004  1=1, Kh 
J  =  J  -  1 
LL  =  k*  -  L  ♦  1 
ISLOClJ)  = 

IF  (K.LT.O)  KTSGH(J)  =  LL 
H I GC  ( J  1  =  R IGC  (  I J 
IF  (f]AK-1)/2)  4<i44, 4090, 4088 
IF  (IAK.Fa.3l  IH1  =  I B 1  ♦  BIGC(J) 

IF  (IAK.F0.4)  Th2  =  T8?  ♦  RlGC(J) 

IF  (K)  SGR  (  J  J  .  E.N.w-BAR  )  B  T  GC  (  J  1  =  C.5*RIGC(J) 

GO  TO  4 oq2 

IF  UkK.fO.l)  TGI  =  TGI  ♦  RTGC(J) 

IF  (IflK.fg.21  TG2  =  TG?  ♦  BIC-C(J) 

IF  (*  T  SGR  (  J)  .F.G.nUAW  )  PIGC(J)  =  P.5*6TGC(J) 

IF  (Kl  SGBfJl  .E'-.l  »  BIGC(J)  =  o.5*6TGC(J> 

COnT  I  mi 
CONTI  1  i  U F 
NTSC  =  NTSC* 

IF  (NDERV.fcU.l.0^.KTYPF.Er.3T  GO  TU  4100 
f)G  =  H*R**3/ld. 

D8  =  T*G 
GC  TO  a  11 0 
OG  -  r«ii*a*«3 
06  = 

IF  (USG.tN.O)  GO  Til  4120 
On  ans  i  =  i, nsg 

OG  =  OG  ♦  A**3*01  IG(  1  )*HMKSGf  II  )/(6.0*FLOAT(NB-1  )  ) 

IF  (nSR.EO.'O  GO  CO  0135 
OU  4  130  I  =  l,wSrs 
L'HX  =  0. 

PCi  4125  .1  =  1,  NRT 

ORX  =  0  4  x  ♦  0220(1, ( J 

v> b X  =  DR x /FL 0 4  T  ( i'J»i  T  ) 

DR  =  !'  R  *•  !)nX*«**S*HJ(KSH(I))/(6.  (>*  FLOAT  (NG-1)) 

COOT  If  !'F 

TGI  =  r  G  1  *  x  L  *  *•  /  (  'G*pLO  A  T  (  Oh4«>- 1  }  ) 
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r 


Tea  Z  TB?.xL.*/tOG.FLOAT(NB*R-i))  Reproduced  From 

TB1  =  tbi*a/(D6*<j*floaT  imbar-i ) )  Best  Available  Copy 

T82  =  T82*A/(0B*<J*FL0AT(MBAR-1) ) 

IF(fcSYMG.EG.O)  TG2=TG1 
IF(NSYMB.EG.O)  T92=TB1 

WRITE (6r4600)  TGl,T8t,  TG? , T62 . DG . OB , XLP4 , XLP5 
«600  FORMAT (lx,8E15.b) 

4700  CALL  PLANE 
CALL  BOLT 

4BO0  K  =  0 

OA=l .o/A 

00  520  0  Irl.r.GT 

00  5200  1=1, 08 I 

IF (NUSF( J. I) .EG.O)  GO  TO  5200 

K  =  K-f  1 

OwOCKJzO.O 
0i»G(K)=0.0 
0*9 (K) =0.0 

IF (NSMEAF.EO.O)  Go  TO  5000 
XSXZ(K)=0. 

XST2(K)=0. 

PESXZG(K1=0. 

PESTZbfM=0. 

FES X2B(*)=0. 

PESTZnCM=0. 

5000  CONTINUE 

DO  5100  M=1,MG 
kf=(«-1  )*f.'GT  ♦  1 
OO  5ir.o  5=1, *8 

IF  ( MUSE  ( f. ,  r*)  •  EG.  0  J  GO  TO  Sjno 
NN=(\-1)*nHT  ♦  J 
FGMN=FGCn.M)*Oa 

0*0(K)  =  FGMN*FFlCM-l*FP5f  f>N)  ♦  pi-.Q  (K  } 

D”*G  (  K  J  =  FGFn*FP2  (  vm  i  *FP5  ('■•r-i  ♦  r> t'  G  f  k  l 
Dr-BfO  =  FGfk*FF  1  (*«*«]  *FPb  i  Nfj )  ♦  CaR(K) 

5100  CONTINUE 
5200  COMI'  OF 
NGNbT  =  K 
L  M  A  X  =  LF6K*NG'JbI 
NGNB  =  0  G  *  N  B 
CUP  =  mplT 

IF  Cr.SwEAP.fc*;.'))  GO  TO  53"0 

XLP8  s  -a.n/ (g.o*FLOAT(«6AP-1  }  *FLOA  T  C.'BAR-l)  *L  } 

XLP7  s  XLPrt/GG 
XLPB  =  XLPF/GH 
SET  UP  “OPt  SM(.Pt5, 

K  =  0 

00  523"  T=l,LR«»AK 
M  =  LP  G  (  T  1 

KM  =  i' 

0  0  5230  J  = 1 , y h A p 
K  =  K  ♦  1 
X  =  G  A  w  (  J  )  *  X  >' 

cx  =  cr-six) 
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K  =  0 

00  52o0  I=1,LSHAR 

h  =  U'B(i)  Reproduced  From 

x  n  r  r.  Best  Available  Copy 

DU  5260  J21.N0AR 
K  8  K  ♦  1 
X  =  BFTR(J)*XN 

cx  =  ros(x) 

SNX  =  SIN (X) 

F  X  73  (*  )  =  Sf.x 
F  X Z a ( K )  =  XN*Cx 
FTZ3CK)=XN*CX 
5260  FTZ«CK)=  -Sf'X*xw**2 
5300  CONTINUE 
RETURN 

R300  FORMAT  f 1H1 ,20X,  33hO  E  P  K  n  P  f  B X  TENDED  N0V„.  1980)/ 

1  I5H«)PANtl  ANALYZED) 

9U0C  FORM#  T ( 7h  FLAT) 

Q500  FORMAT (oh  CURVED) 

ObOO  F URM A  T ( 2?h  metal.  SINGLE  LAYER) 

0700  FURMAT(?aH  PLASTIC.  SINGLE  LAYER) 

OdOO  FORMAT (19b  mETAl.  HONEYCOMB) 

9820  FORMAT  ( 2 1 h  PLASTIC.  HONEYCOMB ) 

Oft/40  FORMAT  PLASTIC,  MOI  T  I  LaYFR) 

9500  FORMAT  (T7n  CLAPPED  -  CLAPPED.  GAMMA  DIRECTION) 

992n  FORMAT  (35"  SIMPLE  -  SIMPLE.  GA**ma  DIRECTION) 

9030  FORMAT  C  3 1  h  FKEfc  -  FREE,  GAU”A  DIRECTION) 

9900  FORMAT  ( 3 6h  CLAPPED  -  SIMPLE,  GAV.m*  DIRECTION) 

9050  FURMAT  (39H  CLAPPED  -  FREfc.  GAMMA  OTRECT ION) 

0055  FORMAT (31 h  Tup  SPRINGS.  Gamma  DIRECTION  ) 

0060  FORMAT  (3bH  CLA*-PEi)  -  CLAf-PFD,  BETA  DIRECTION) 

9980  FORMAT  (3/ih  SIMPLE  -  SIMPLE.  BETA  DIRECTION) 

9990  FORMAT  (36m  FREE  -  FREE,  BETA  DIRECTION) 

0000  FORMAT  (35h  CLAPPED  -  SIMPLE,  BETA  DIRECTION) 

0050  FORMAT  (33b  CLAmPE.0  -  FPEF  ,  BETA  OIREC  flON) 

f) . » 5  S  FORMAT  (  3  o  h  T.jh  SPRINGS.  BETa  DIRECTION  ) 

0100  FURMAT  (26M0PESPtiNSh  OPTION  -  RI.ASTIC) 

0200  FORMAT  (  3aH0RESRU:''SE  OPTION  -  ELASTIC-PLASTIC) 

0250  FORMAT  (  ?  6  R*  C  I  n  R"  R  !  I  a  COUPLING  INCLUDED) 

OdO0  FORMAT f ! 7HCS  IPuC ! OKAL  MODEL/ 

1  a?H  NUMbER  uF  ga^ma  MQf  ES  (MG)  =  13/ 

P  97b  mi-'EER  ub  ^EIA  MOOES  f'-8)  =  13/ 

3  a  7  h  0  NUMbE  h  iJL  G  A  A  INTEGRATION  HO  T u T  S  (HOAR)  z  13/ 

a  U7h  HF>  UF  b E  T A  INTEGRATION  POI-jTS  (N8«.R)  r  13/ 
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Reproduced  From 
Best  Available  Copy 

5  «7H  NUMBER  i)F  Z  INTEGRATION  POINTS  CL6AK)  s  13) 

0820  FORMAT  (24HOMOOAL  COMBINATIONS  USED) 

0830  FORMAT  (3X,2I«) 

0850  FORMAT  (35H0  LENGTH  OF  PANEL.  IN  (XLP)  =  Elb.8) 

0900  F0RMAT(35H  *IDTH  OF  PANEL,  IN  (THETaO)  =  Elb.8) 

1000  FORMAT ( 35H  SU8TEN0E0  ANGLE.  DEG  (THETAO)  =  Elb.8/ 

2  35H  RADIUS,  IN  (A)  =  E16.6) 

1100  FORMAT ( 35H  THICKNESS.  IM  =  Elb.8/ 

1  35H  DENSITY,  i.8-SEC**?/lN**a  =  Elb.8/ 

2  35H  ELASTIC  MODULUS.  P3I  =  Etb.8/ 

3  35H  POISSON  IS  RATIO  =  Elb.8/ 

a  35H  YTFLO  STRESS,  PSI  =  Etb.8/ 

5  35H  STRAIN  nARuENlKG  SLOPE,  PSI  =  EIb.8/ 

b  35H  ULTIMATE  strain,  I*. /In  (EPSIF)r  Elb.8) 

1200  FORMAT ( 2 b H o I n  1  T 1 AL  I MPERFt  CT IONS.  If*/ (5E14.6) ) 

1  Juft  FORMAT  f 1 7H0 T IMt  I  ■  -  E  U  R  M  A  T 1 0  N / 

1  «2H  INTEGRATION  STEP  SIZE.  SEC  (OElTIm)  =  Elb.8/ 

2  a2H  STOP  TI-E,  SEC  (TSTOP)  =  EIh.h/ 

3  «£H  PRINT  FRt'J'iENCY  (PRINT)  =  Elb.8) 

2050  FORMA  I  l<i o HOC OCRUl NATE  SURFACE  POSITION  (HBAR) ,  I  u  Etb.8/ 

1  1  3m 0 LATER  nin-hEk, 2?X,al  IS/ (  31  x.al  15) ) 

2100  FORMAT  (27h  CUMULATIVE  THICK  f- ESS,  T.-> ,  1  3  X  ,  bE  1  5  .  h  T 
2200  FORMAT  (3?h  '-aSS  OEnSJTi,  L  B-SEC  *  *2/ IN  **«  ,  8X  ,  fcE  )  5 .  b  T 

23>)0  FORMAT  (  33h  -UUUL'JS  OF  ELASTICITY  -  X,  PS  I  ,  7  X  ,  bE  1  5 .  b  ) 

2«00  FORMAT  (aOH  *-uf)ULUS  OF  ELASTICITY  -  THETA,  PST  tblS.b) 

2500  FORMAT  {??m  PU1SS0«(S  RATIO  -  X . 1 d * . &E 1 5 . b  ) 

PbOO  FQRmat  (?^jH  PUISS'WS  RaTIO  -  T  mf  T  A  ,  1 «  X  ,  bE  1  5 . 6  ) 

2700  FORMAT  (31h  TEwbILE  ULTIMATE  STRESS,  PSI ,  9X ,  HE  1  5  .Vj 

2800  FORMAT  (35*  COMPRESSIVE  L'LTJ-ATt  STRESS,  PS  J  ,  5*  ,  bt  1  5 .  b  ) 

2900  FORMAT  (2*h  (L'OlLE  YlfcLP  STRESS,  PS  1 ,  1 2  X  ,  bE  1  5  .  t, ) 

3000  Format  orh  COMPRESSIVE  YIELD  STRESS,  PS  I  ,  8X  ,  bfc  1  5 . 6 ; 

2850  FORMAT  (25 M  i-tAH  MODULUS,  PS I , 1 9X , bt 1 5 . b ) 

310  0  FORMAT  (bPHOCnPt  MODULUS  CF  ELASTICITY  PARALLEL  TO  CORE  DEPTH  (ECi 
1,  PSI  =  E15.K/ 

2  b 3h  Shear  modulus  UF  CDR'F  (GO,  PSI  = 
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Reproduced  From 
Best  Available  Copy 

3  £14.6/ 

4  62H  CORE  CELL  SIZE  (DC).  1 N 

5  E15.6) 

3200  FORMAT  (41H0**  rtAKNINS  «*  INCONSISTENCY  IN  SYMMETRY} 

3300  FORMAT  (4?H0  CUOROI'IATE  SURFACE  POSITION  (H0AR),  Iw  El6.fi/) 

3400  FORMAT  (3AH0NU*iflEH  QF  DISCRETE  ELASTIC  SPRINGS  “  13/ 

1  52HONIJWBER  OF  6UIJN0ARY  POINTS  *ITH  TORSIONAL  SPRINGS  =  13) 

3500  FORMAT  (37HO  Iu-PLANE  TU-V)  EDGE  CONSTRAINTS  -  ) 

3b 00  FORMAT  (21H  CLAMPED  -  CLAMPED,) 

3700  FORMAT  (15H  FREE  -  FREE,) 

3800  FORMAT  (18H  CLAmpeo  -  FREE.) 

3900  Format  (1 H* , ?9X , \ 5mG AMMA  DIRECTION) 

aOoO  format  ( 1 h^,?ua, 1«mHETA  DIRECTION/ 

1  33H0  OtiT-OF-PLA^e  (A)  CONSTRAINTS  -  ) 

41. )0  FORMAT  (3bM0SHEAK  DEFORMATION  OF  CORE  NEGLECTED) 

<420°  FORMAT  (3SHCSHEAR  D€  F  liRMA  T  J  On  OF  CORE  INCLUDED) 

END 

•DECK  S I GM a 

SU8R0NT  I  Nf  SIG.J«  (I.J.m) 

This  y/o  0ETERh1<ES  THE  STRESS-STR  A  Ik'  RELATIONSHIPS  FOR 
ELASTIC  A  NO  /  (JR  PLA3TTC  RESPONSE. 

Subroutine  completely  revised  mapch,  1976. 

K  -  JNDFX  'JF  THE  I'lTEQPAT  TDtv  POINT  Jm  ThE  z  DIRECTION. 

I  *•  I'JDtX  *jF  Tit  INTEGRATION  POT!,T  IN  THE  GET  a  DIRECTION. 

J  -  INDEX  OF  The  INTEGRATION  PUT?  T  in  The  gamma  OIhECTIUN. 

★CALL  C4LK1 

•Call  C8LK3 

•  CALL  C9LK«J 

•  CAUL  CPl.Kh 
•CALL  CSLk7 

•  CALL  Cr.OVA 
•CALL  CBLANk 
•CALL  CPLK15 


DATA  TOL/5.0E-3/ 

IF (IF ! RST.GT.O)  GO  TO  300 
IFIFST  =  1 
EP0  =  s  tgo/el 
t  p  P  =  n  .  0 
CNR  =  0.'! 

SIGN?  =  3  f  G  *J  *  •  d 


I2S 


CNl  =  (0.5  -  TNUJ/EL 

CN3  =  1.0/EL 

TNUSQ  =  TMU**2 

CNfl  =  (1,0  -  TNUSUJ**2 

CN12  s  (1.0  -  TrtU  <►  TNUSQ}  /CN4 

CN13  =  U.O  -  4.0*TtJU  ♦  TNUS'))/CQ4 

CN4  =  0.75/( (1.0  ♦  INU)**2) 

CN6  s  EL/fl.O  «  TN'J**2) 

CN7  =  EL#0 .5/  ( 1 r  0  ♦  TNU) 

CN5  =  1.0/CN7 
LC  *  0 
LCMAX  a  100 
00  100  L=t,L«AX 
IF  (NELP.E0.1J  GO  TO  70 
ALXX(L)  =  0.0 
ALTT(L)  =  O.o 
ALXT(L)  a  O.o 
BEl(L)  =  0.0 
8F2 (L 1  =  0.0 
8E3CL)  =  0.0 
TTNU(L)  a  Tnu 
70  EPBU(L)=0.0 

100  KT(L)  a  l 

300  K  SlJM  a  0 

IJ  =  LBA9*(M-1J 
00  3050  Kal.LBArf 
L  a  IJ  ♦  K 
HI  =  ZFC<) 

0ETEP:«IME  appropriate  REGION. 

KEY  a  «r(L) 

IF(kE*.GT.3)  GO  TO  350 
GO  TO  (4 no ,b0o, / 0 0 J ,  KEY 
350  IF  (  (KEYM  ) /a.E-J.Kt  Y/2)  GO  T<)  609 
GO  T O  700 

*<FGTU-J  1.  ELASTIC  CURVE. 

400  K  SlJM  Z  K  S  t  M  ♦  l 

IF  (KS‘JM.GT.1  )  GO  TO  450 
01  a  C  N  6  *  (FXX  ♦  r,^ij*ETT) 

02  s  CNf?*f£TT  t  T > , u * E X X ) 

03  a  CN7*6XT 

04  a  CN6*fXKXX  ♦  TwUaXKTT) 

05  =  CNb*f’XKTT  f  rNU*XKXX ) 

06  a  CM7*XKXT 

IFfNHlJCK  .EU.  9  .UR.  iVC  ALL  .GT.  0)  GO 
IF (02  .GT.  BUCK)  GO  TO  «30 
02  a  BUCK* (8UCK/U2) .ft 
430  CONTMiJk 


S1A(P) 

— 

01 

S2A(M> 

s 

D2 

S3A (w) 

zz 

^3 

S4A  f v) 

r 
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S5A (m j 

= 

n  5 

Reproduced  From 
Best  Available  Copy 


TO  430 
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62  s  02  ♦  H 1 *05 

63  s  03  ♦  HI *06 


Reproduced  From 
Best  Available  Copy 


SIGBO  =  G 1  * (G 1  -  G2 )  ♦  G2**2  +  3.0*G3**2 
IF  (NELP.E0.2)  GO  TO  470 
EP80(L)=SIG80 
GO  TO  3000 

470  IF (SIG8D.GE.3IG02)  GO  TO  500 
EPBO(L)  a  S IG80 

GO  TO  3000 


50n  KT(L)  =  KEY  ♦  l 


LINEARLY  INTERPOLATE  on  sigma  bar  to  correct  for  overshoot. 
SUSIG  =  SORT (SIGBO) 

B2  =  SORT ( E  P  B  0 ( L  J ) 

B 1  =  CSIGO  -  62) / (3QSIG  -  R2) 

G 1  =  SX  X  CL)  ♦  H 1  * (G  1  -  S X  X  (L )  ) 

G2  s  STTCL)  ♦  6 1  * ( G2  -  STT  CL ) ) 

G3  =  SXTCL)  +  8 1  * ( G3  -  SXT(L)) 

SIGXXl CL)  =  G1 

SIGTTICL)  =  G2 


SIGXTt  CL)  =  G 3 

T  1  =  CM  ?*  (G1  -  T i'J U * G 2 ) 

T2  =  CN3*CG2  -  )Nu*Gl) 

T  3  =  CN5*G3 

EXX1CL)  =  T1 
ETT1CL)  =  T  2 
£  X  T  1  C  L  )  s  T  3 

EPBO  =  SGRHC<',12*(ri**2  ♦  T?**2)  -  CN13*T1*T2  ♦  CN0*T3**2) 

EPBCL)  =  EPBO 

EPBO CL)  =  EP60 

IF  (JFIRST.EU.OJ  JFIRST  =  1 

GO  10  300O 

REGIOhS  2  AMO  4.  PLASTIC  LOAOtuG. 

bOO  E P 9 D P  =  EPBCL) 

H2  =  EXX  *  H  J  *  xK  a  X  -  R  E  1  (  L ) 

H3  =  ETT  +  HI  *  X  K  T  T  -  BE?  CL ) 

hu  =  EXT  +  hiaX'XaT  -  0E3CD 

CM?  =  TTnijCL) 

II  =  0 

blO  II  =  II  +  1 

CM2?  =  C^?**? 

EPBO  =  SUWTCCCl.O  -  cm?  ♦  c  ••*??)*  C  02**2  *  H3**2)  - 
l  Cl. 5  -  U.O*Cw?  *  CM??)*H?*HS)/Cl .0  -  CM??)**?  + 

?  0.7S*hu**2  /  ( 1  *i'  +  CM?)**?) 

OELEP  =  F PHD  -  EPHOIL) 

EPP  =  (FP*OFLEP  ♦  EL  *EPBO (L ) ) /EPBO 
IF  (Till). GT. 0.0 )  CM?  =  .5  -  EPP*CM  1 
IF  CA3SCCM?-TTMtJlL)  )  .LT.n.OonS)  GO  TO  b20 
IF  (Il.GT.20)  GO  TO  olS 
TTMI'CD  =  Cm? 

GO  TO  bl»> 

blS  ‘'•RITE  Ch.SSOO)  C !  2  »  f  r  M 1 1  f  L )  .  T  1  ’  F. 

GO  TO  -l  1  ni) 
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620  CN2  =  TTNU(L) 

IF  (EP80.LE.EP60P)  GO  TO  6S0 
630  EPS (L)  =  EPBO 

IF  (EPP.GT.EL.OR.EPP.LT.EP)  go  to  ao 00 

51  =  EPP/n.O  -  CN2**2) 

52  s  0.5*£PP/(1.0  ♦  CN2) 

G1  =  Sl*fH2  ♦  CM2*H3)  4-  A  L  X  X  f  L  ) 

G2  =  St  *  CH3  ♦  C  n  2  *  H  2 )  ♦  ALTT(L) 

G3  =  S2*H«  ♦  ALXT (L) 

GO  TO  3000 

SECOND  TEST  FOR  UNLOADING  IN  EITHER  REGION  2  OR  0  . 
650  Ol  =  Xl  A  (M)  -EXX1  (L)  ♦  H1*X4A{M) 

Q2=X2A (M).ETTl(L)  +  Ht*X5A(M) 

I)3  =  X3AfM)-tXTl  (L)  ♦  Hl*X6AfM) 

IF  CEP* ED, 0.0) GO  TO  660 
PUSXXCD-SIGXXUL)  +  ALXX(L) 

P2=STT  (L)-SIGTTHL)  +  ALTTIU 
P3=SXT(D-SIGXTUL)  +  AL  X  T  (L  ) 

GO  TO  670 
660  P1=0.0 
P2  =  0.0 
P3=0.0 

610  E 1 s 0 1  -  M2 
E2  =  D2  *  m3 
E3=U3  m  hu 

Gl  =  Pt-CN6* (Et  +  r<u*£2) 

G2sP2-CN6*(E2+rNU*El ) 

G3=P3-CM7*E3 

MsGt-Pl 

A2=G2-P2 

A3=G3-P3 

SIGHO  =  Al*fAl-A2)-»-A2**2  +  3.0*A3**2 

IF  (STOOD. GE.SI302. AND. DELEP.GE. 0.0)  GO  TO  630 

KYCL)=KEY+1 
TTNU(I.)sT-JU 
£PBO(L )=SIGBD 
dE  1  (L)=Q1  +  *it  1  ( L  ) 

Q  E  2  ( L  )  *  Q  2  ♦  3E2CL) 

BE3(L)=U3  ♦  H  t  3  ( L  ) 

IF (EP.EQ.O. )G0  TU  3000 
ALXX (l)=P1 
ALTT (l)=P2 
ALXT (L)=P  3 
GO  TO  3000 


Reproduced  From 
Best  Available  Copy 


REG  I DO  3.  ELASTIC  UNLOADING  -  RELOADING. 


700  El 

z 

3 E 1  CL)  -  EX x 

E2 

z 

BE 2 ( L  )  -  ETT 

E  3 

2 

H E 3 ( L )  -  EXT 

ct 

s 

ALXX  fL) 

C2 

z 

ALTTfL) 

C  3 

= 

ALXT  fL) 

G  1 

2 

r.  1  - 

H  l  ★  X  K  X  X 
H1*VK  TT 
H| *  X  K  X  T 


Cn6*(E1  ♦  T 1,1 1 1  *  1 2  ) 
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CN6*(E2  +  TWINE  1 ) 
CN7*£3 


Reproduced  From 
Best  Available  Copy 


G2  =  C2  - 

S3  s  C3  ■  LiMf*t3 

M  =  SI  -  Cl 
A2  s  G2  -  C2 
A3  =  S3  -  C 3 

SIG8D  =  A 1  * ( A 1  -  A2)  ♦  A 2**2  ♦  i.0*A3**2 
IF(SIGBD,GT.3IG02)GO  TO  800 
EPB0CL)=SIG80 
GO  TO  3000 


LINEARLY  INTERPOLATE  ON  SIGMA  BAR  TO  CORRECT  FOR  OVERSHOOT. 

800  82  s  SORT (EP60 (L) ) 

SUSIG  s  SORT (SIGbO) 

IF  t82.GT.SlG0)  GO  TO  800 
NC  =  U 

820  81  =  CSQSTG  -  SIGU)/(S«$IG  -  P2) 

NC  =  NC  +  1 
IF (NC.GT.5)  GO  TO  830 
DEL  1  =  R 1  *  t G 1  -  SXX (L) ) 

0EL2  s  HI « (G2  -  SI  T (L)) 

OEL3  =  3 1  *  (G3  -  S<T(D) 

G 1  =  G 1  -  0EL1 
G2  =  G2  -  DEL 2 
G3  =  G3  -  0EL3 
A 1  -  G 1  •  ALXX(L) 

A2  =  G2  -  ALrrtL) 

A3  =  G3  -  ALXTCU 

SOS I G  s  SQRT(A1*(A1-A2)  +  A2**2  ♦  3. 0*A3**2) 

IF (A83(SQSIG-SIG0) /S1G0.GT.  IOL)  GO  TO  820 

GO  TO  835  t 

330  WRITE  (fe,  57 00  )  NC  #  K , I , J , K  E  Y , SO  $ I G , 8 l , 62 »  T I 
LC  =  LC  ♦  1 

IF  (LC.GT.LCMAX)  GO  TO  0100 
835  CONTINUE 

0EL1  =  Gl  -  Sxx(L) 

DEL2  =  G2  -  ST  I  CL) 

DEL  3  =  G  3  -  SxrtU 

T 1  =  XlA(w)  +  CW3*(0EL1  ”  T  N u  *  D  E  L  2 )  ♦  H1*X0A(mj 
12  s  X  P  A  (  m  )  +  C  i*i  3  *  ( 0  E  L  2  -  T-\H*nELl)  ♦  H1*X5A(M) 

T3  s  X 3 A ( m )  +  CNb*0£L3  +  nl*XbAH) 

GO  TO  880 

San  WRITE  (6,5200)  K ,  I ,  J  ,  KEY  ,  TIf-E  .  R?t  SOS  (G 

T 1  =  Exx  +  H  1  t»  X  *  X  A 

T 2  =  £ T  T  +  HI*  AMT 
T3  =  EXT  +  H  l  *  A  i\  X  I 
LC  =  LC  ♦  1 

IF  (LC.GT.LCMAX)  Go  TO  a  100 
380  EXX1 (L)  =  Tl 
ETT1  (|.)  =  12 
EXTl(L)  =  F3 
H2  =  ft  -  8EI (L) 

H3  =  T2  -  8E2(L) 

Ha  =  T3  -  mE  3 ( L j 

fcPRO  =  S0WT(CU12*  IH2**2  +  H?**2)  -  Cf'- 1  3*H2*H  3  + 

£PRO(L)  =  FHiiH 
EPHlL)  =  EPaD 
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SIGXXHL)  =  G l  Reproduced  From 

SIGTTl(L)  =  G2  Best  Available  Copy 

SIGXTHU  =  G3 

K  Y  (L)  =  KEY  +  1 
GO  TO  3000 

3000  SXX(L)  =  G 1 

STT(L)  ~  G2 

SXTCL)  =  G3 

3050  CONTINUE 

KSUMA(M)  s  KS'JM 

RETURN 

ERROR  RETURN. 

4000  WRITE  (fa,  5300  )  EPP .  K  ,  I .  J ,  T  IMF ,  £PBf) ,  EP80P »  EP30  (L ) 

4100  WRITE  (6.5400) 

KERR  s  1 
RETURN 

5200  FORMAT  (22H  IMMEDIATE  RELOADING  ,4I3,3E15.fa) 

5300  FORMAT  (23H0EPP  13  OUT  OF  RANGE,  EPP  =  E14.6/ 
l  3I3.4E15.6) 

5400  FORMAT  (21  HO  SOLUT I  UN  IS  UNSTABLE) 

5500  FORMAT  (2feH  VALUE  OF  NtJ  ;n0NT  CON  VERGE , 2E 1 5 . fa . 1 5H  TImE,  SEC 

1  £15.6) 

5700  FORMAT  ( 3 AH  CAN  nOT  TOTALLY  CORRECT  FOR  OVERSHOOT /S  15, 4E  15. 6) 
END 

*OECK  SIGMA9 

SUBROUTINE  SIGMAH  (K 1 ,K2, J, I ,f») 

THIS  PROGRAM  COMPUTES  INELASTIC  STRESSES  FUR  STIFFEwERS. 

K 1  -  gamma  STIFFENER. 

K2  -  BETA  STIFFENER. 

J  -  BETA  POINT  HUMBER. 

I  -  gamma  POInT  UijmBER. 
m  -  GRID  point  number. 

ONLY  ONE  STIFFENER  IS  ANALYZED  PER  CALL. 

•CALL  CBLK1 
•CALL  CBLK7 
•CALL  CBLK15 
•CALL  C8LK16 
•  CALL  CHLAfK 

IF  (NF IRST.GT.O)  GO  TO  300 
NF IRS  T  =  l 

MAX  =  NSMAX* (NSGmB  ♦  USBaMB ah) 

DO  100  Lsl.MAX 
ALX(L)  =  0. 

BE  X (L )  =  0. 

100  KYX(L)  =  1 

500  K SUM  =  0 

IF  (Kl.EO.O)  GO  TO  320 
GAMMA  DIRECTION. 

IJ  ~  (K1— 1)*^BAK  ♦  I 
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Reproduced  From 
Best  Available  Copy 


NS  s  MSEGG(Kl) 

KKK  S  K 1 
ELX  =  ESTRGCK1) 

SIGX  s  S IGOGT  C K  1 ) 

EPOX  =  fcPGG(Kt) 

EPX  s  ETSTRG(Ki) 

ex  =  exx 

XK  =  XK  XX 
Xic  =  XIAfM) 

X4C  =  X4A(«) 

GO  TO  330 
beta  direction. 

320  IJ  s  NSGMB  ♦  (K2-t )*\BAR  ♦  J 
NS  =  NSEGD (K2) 

KKK=NSG* (K?-l ) aYBAR+J 
ELX  =  ESTR6(K2J 
SIGX  =  S  f  GOBT ( K  2 ) 

EPOX  s  FP0d(*\2) 

EPX  =  ETSTRB(K2J 
EX  =  ETT 
XK  =  XK  T  T 
XIC  =  X?A(^) 

X4C  =  X  5  A ( M ) 

330  LO  =  NS MAX  *  (IJ-l) 

DO  3030  Ksl.NS 
L  =  LO  +  K 
HI  S  /FHCK,KKK) 

KEY  5  K  Y  X  f L ) 

IF  (KEY. EG. 1)  GO  TO  400 
IF  ( (KEY-M  ) /2.E0. KEY/2)  GO  TO  600  i 

GO  TO  700 

REGION  l.  ELASTIC  REGIO-J. 

400  K SUM  =  KSlJM  +  1 

IF  (K3UM.GT.1)  GU  TO  450 
01  =  E  X  *  F  L  X 
04  =  < K  # E L  X 

450  G 1  =  01  ♦  M 1  * 0 4 

IF  (NELP.E0.1I  GO  TO  3000 
IF  (ArtS(Gt)  .l-T.oIGX)  GO  TO  3000 
INITIAL  YIELDING. 

K  Y  x ( L )  =  KE  Y  *  1 
G2  =  SIGX 

IF  (Gl.LT.O.O)  62  =  -SIGX 
G 1  =  G? 

SIX  1(1.)  =  61 
E  X  1  ( l. )  =  G 1  /ELX 
IF  (MF IKST.EN.OJ  mFTRST  =  1 
GO  TO  3000 

REGIONS  2  AND  4.  INELASTIC  REGIONS.  > 

600  HR  z  EX  +  HI  *A'\  -  4£  X  (l_ ) 

H2U  =  X  1 C  +  Hl*x4C  -  BE  X  (I. ) 

H  3  =  AHS(H2) 
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IF  (H3.LT.ABS(ri2U)  )  GO  TO  fe50 

OELEp  s  H3  -  £POX 

EPP  s  (EPX*D£LEP  ♦  3IGXJ/H3 

61  =  EPP*H2  +  ALXCLJ 

GO  TO  3000 

UNLOAD. 

650  8£X(L)  =  2.0*8£XIL)  ♦  H20  -  E  X  1  (L 1 
1  ALX(L)  s  ALX (L  i  ♦  SX(L)  -  SIXi(L) 

G 1  =  (EX  ♦  H1*XK  -  BEX  (L) ) *&LX  ♦  ALX (L) 
KYX(L)  5  KEY  t  1 
GO  TO  3000 

REGION  3.  ELASTIC. 

700  H2  s  E X  ♦  Hl*XK 

A 1  s  (H2  -  6EX(L))*£LX 
G t  a  A 1  ♦  ALX(U 

IF  (AHSC At ) .LE.SIGX)  GO  TU  3000 
YIELD.  % 

SXL  =  SX(L) 

81  «  ABSfSXL  -  ALX (L ) ) 

ea  =  ( s i g x  -  an/uasui)  -  01) 

G 1  =  SXL  +  b2« (G 1  -  SXL) 

SI XI  (I.)  =  G  1 

E  X 1  ( I. )  =  X1C  ♦  HJ*X4C  ♦  (61  -  SXL)/ELX 
KYX(L)  =  KtY  +  1 

V 

5000  SX(L)  =  Gt 
5050  CONTINUE 

t  KSIJM8  ( I J  5  =  K.SL'M 

3100  RETURM 
END 
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